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INTRODUCTION  AND  SUMMARY 

\ 

This  report  is  concerned  with  two  nonlinear  propagation  effects 
which  limit  the  transmission  of  intense  laser  radiation  through  the 
atmosphere.  One  effect  is  that  of  self-defocusing  due  to  heating, 
which  can  cause  an  undesirable  de  rease  in  the  delivered  power  density. 
The  defocusing  results  from  photon  absorption  leading  to  atmospheric 
heating  and  consequently  to  a  local  decrease  in  the  refractive  index; 
the  result  is  a  lens  effect  which  causes  the  beam  to  diverge.  This 
instability  is  apparently  significant  even  for  exceedingly  low  ab¬ 
sorption  coefficients.  Consequently  it  is  necessary  to  have  reliable 
information  in  the  far  wings  of  atmospheric  absorption  lines  where 
laboratory  data  is  meager,  and  where  conventional  pressure  broadening 
theories  are  inapplicable.  This  report,  then,  presents  theoretical 
studies  on  certain  aspects  of  very-far-wing  collision  broadening, 
which  are  particularly  pertinent  to  the  problem  of  thermal  defocusing 
in  the  atmosphere. 

The  second  phenomenon  studied  is  the  stimulated  Raman  effect. 

This  work  represents  a  continuation  of  previous  research'*'  which  was 
concerned  with  catastrophic  beam  depletion  resulting  from  stimulated 
Raman  scattering  in  the  atmosphere.  The  new  work  deals  with  appli¬ 
cations  of  the  previously  derived  general  results  to  an  analysis  of 
the  generation  and  amplification  of  short  intense  pulses  of  Raman 
backscattered  light. 

Chapter  1  describes  some  preliminary  theoretical  studies  on  the 
far-wing  collision  broadening  of  the  rotational  spectrum  of  a  polar 
molecule.  The  studies  were  motivated  by  the  currently  prevalent  belief 
that  a  substantial  part  of  the  continuous  background- absorption  by  the 
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atmosphere  in  the  wavelength-vicinity  of  lOy  is  due  to  the  very-far-wing 
collision  broadening  of  the  rotational  spectrum  of  ^0. 

The  basic  approach  of  the  herein  described  studies  may  be  regarded 
as  an  extension  of  the  classical  Debye  theory  of  dielectric  relaxation. 
In  this  theory,  the  origin  of  dielectric  absorptivity  arises  from  the 
random  interruption  of  the  rotational  motion  of  the  individual  molecular 
dipoles  by  impact-type  collisions,  i.e.,  collisions  of  infinitesimal 
duration  —  with  other  atoms  or  molecules  of  the  medium.  In  the  present 
work,  the  theory  has  been  extended  to  take  account  of 

(a)  the  fact  that  (in  contrast  to  Debye's  original  as¬ 
sumption)  the  frequency  of  the  incident  light  wave 
is  large  compared  to  rotational  frequencies, 

(b)  quantum  corrections  arising  from  the  fact  that  -fiu)  >  kT, 
and, 

(c)  the  actual  finiteness  of  time  duration  of  collisions. 

While  this  work  was  nearing  completion,  the  authors  became  aware 
of  very  recent  data  of  D.  Burch  on  laboratory  studies  of  ^0  - 
mixtures  in  the  lOy  region.  These  studies  indicated  (cf.  Summary  of 
May  16,  1967  meeting  on  Infrared  Atmospheric  Absorption,  page  2)  that, 
despite  the  fact  that  the  relative  concentration  of  1^0  in  the  studied 
mixtures  was  <  2%,  the  collisions  responsible  for  absorption  were  those 
of  ^0  molecules  with  each  other!  Such  a  result  was  found  to  be  un¬ 
explainable  within  the  framework  of  the  above-described  theoretical 
studies;  moreover,  the  magnitude  of  the  experimental  absorption  was 
found  to  be  several  orders  larger  than  the  theoretical  prediction. 

The  first  attempt  to  remove  the  discrepancy  consisted  in  taking 
account  of  the  shortening  of  the  time-duration  of  collisions  via  the 
mechanism  of  attractive  forces  (associated  with  dipole-dipole  inter¬ 
actions);  as  shown  by  the  detailed  expressions,  such  shortening  would 
lead  to  an  enhancement  of  the  far-wing  absorption.  However,  with 
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physically  reasonable  values  of  the  relevant  molecular  parameters,  only 
marginal  improvement  was  obtained;  the  calculated  absorption  still  fell 
short  of  the  experimental  value  by  almost  three  orders  of  magnitude 

Under  these  circumstances,  it  was  felt  necessary  to  introduce 
some  essentially  new  physical  ingredient  into  the  theory.  One  such 
possibility  has  suggested  itself  and  is  presently  being  pursued.  Briefly, 
it  is  now  felt  that  the  origin  of  the  far-wing  absorption  in  1^0  vapor 
arises  from  the  angular  dependence  of  the  (strong)  dipole-dipole 
interaction.  Namely,  at  sufficiently  close  —  yet  still  accessible  — 
distances  of  approach  of  two  colliding  ^0  molecules,  this  angular 
dependence  provides  a  potential  trough  for  the  rotational  motion  of 
each  molecule;  within  this  trough,  the  molecules  execute  a  type  of 
hindered  rotation  —  or,  as  it  is  commonly  designated,  a  "librational" 
motion.  The  decisive  point  is  that,  according  to  preliminary  estimates, 
the  frequencies  associated  with  this  librational  motion  are  substantial¬ 
ly  higher  than  those  of  free  rotation  —  and,  in  fact,  may  be  expected  to 
encompass  the  1000  cm  (10p)  region.  Under  these  circumstances,  it 
turns  out  that,  according  to  the  general  theory  of  collision  broadening, 
the  occurrence-probability  of  these  frequencies  —  in  particular,  those 
in  the  immediate  vicinity  of  the  external  frequency  —  becomes  the 
principal  factor  in  determining  the  absorption.  Preliminary  calculations 
indicate  that  the  resultant  absorption  far  outweighs  that  obtained  in 
our  previous  theoretical  calculations. 

In  conclusion,  then,  while  it  is  too  early  to  be  certain,  the 
hindered-rotation  mechanism  looks  presently  rather  promising. 

Chapter  2  contains  a  discussion  of  pulse  amplification  by  means 
of  the  stimulated  Raman  effect  and  of  a  second  related  topic,  namely, 
the  generation  of  intense,  ultra-short  Raman  backscattered  pulses 
recently  observed  in  the  laboratory.  The  first  subject  in  Chapter  2 
is  the  Raman  amplifier  which  consists  of  a  length  of  Raman  active 
material  through  which  both  a  beam  of  laser  light  and  a  beam  of  light 
at  the  first  Stokes  frequency  of  the  laser  light  in  the  material  are  passing. 
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The  beams  are  traveling  in  opposite  directions  and  hence  the  second  beam 
(the  Stokes  shifted  beam)  will  be  amplified  by  the  stimulated  emission 
of  backscattered  light  from  the  first.  For  an  input  Stokes  wave  of 
arbitrary  time  dependence,  the  time  dependence  of  the  output  has  been 
formulated  in  terms  of  the  solution  of  a  differential  difference  equation. 
The  solution  of  this  equation  has  been  reduced  to  quadratures  for  an 
arbitrary  input,  and  has  been  obtained  in  closed  form  for  the  special 
case  of  a  step  function  input. 

The  Raman  pulse  generator  is  a  device  which  was  first  investigated 
by  Maier,  et  al.2  In  their  experiment  a  carbon  disulfide  cell  was  il¬ 
luminated  by  a  beam  of  intense  laser  light.  Observations  at  the  entrance 
of  the  cell  revealed  intense,  short  pulses  of  backscattered  Stokes 
radiation.  They  have  suggested  that  "A  probable  mechanism  for  the 
initiation  of  the  pulse  is  the  abrupt  onset  of  backward  stimulated 
Stokes  emission  near  the  exit  cell  surface,  accompanying  the  occurrence 
of  laser  self-focusing  in  that  region."  As  a  support  for  their  sug¬ 
gestion  we  have  constructed  a  specific  model  of  the  influence  of  self- 
focusing  on  the  amplification  of  spontaneously  emitted  Raman  waves. 

This  model  predicts  the  height  and  shape  of  the  emitted  pulse  in  terms 
of  the  solution  of  a  nonlinear  differential  difference  equation.  The 
equation  is  currently  being  solved  by  computer;  it  is  anticipated  that 
the  analysis  will  be  completed  by  the  end  of  the  present  contract  period. 


2.  M.  Maier,  W.  Kaiser  and  J.A.  Giordmaine,  Phys.  Rev.  Letters,  _x7_, 
1275  (1966), 
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INTRODUCTION 


This  report  describes  some  preliminary  model-type  theoretical 
studies  on  the  far-wing  broadening  of  the  pure  rotational  spectrum  of 
a  polar  molecule.  The  studies  were  motivated  by  awareness  of  the  preva¬ 
lent  belief1  that  a  substantial  part  of  the  continuous  background- 
absorption  in  the  atmosphere  in  the  wavelength  vicinity  of  lOu  is  due 
to  the  very-far-wing  broadening  of  pure  rotation  bands  of  HjO  with 
maximal  intensity  at  wavelengths  v50p.  As  will  be  described  more  fully 
below,  the  basic  approach  is  an  extension  of  that  employed  in  the  cor- 
responding  problem  of  the  far-wing  broadening  of  near  infrared  vibrational 

lines  reported  previously. 

While  this  work  was  in  progress,  the  authors  became  aware  of  some 
very  recent  data  of  D.  Burch  (Philco)3  on  laboratory  studies  of  H.,0 
absorption  in  the  lOp  region,  which  indicated  that  the  basic  broadening 
agency  is  self-broadening,  i.e.,  HjO  -  »2 0  collisions.  Such  a  result 
does  not  find  any  natural  explanation  within  the  framework  of  the  above- 
C  described  theoretical  studies;  moreover,  preliminary  comparisons  indicate 

that  the  absorption  observed  by  Burch  is  almost  three  orders  of  magnitude 
larger  than  could  be  accounted  for  by  these  studies.  A  drastic  modifi¬ 
cation  of  the  theoretical  approach  is  thus  required.  Such  a  modification 
has  in  fact  suggested  itself  to  us,  but  its  development  is  presently  in 

its  very  initial  stage. 

Under  these  circumstances,  we  feel  it  appropriate  to  devote  the 
Ml„  body  of  this  report  to  a  description  of  the  studies  which  have  been 
carried  out,  and  of  their  comparison  with  Burch's  data,  which  reveal 
their  inadequacy.4  The  report  then  concludes  with  a  brief  description 

of  the  new  approach. 


Q 
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1. 


DISCUSSION  OF  THE  MODEL 


The  specific  model  on  which  the  calculations  of  this  report  are 
based  may  be  regarded  as  an  adaptation  of  the  classical  Debye  theory  of 
dielectric  relaxation.  In  this  theory,  the  collision-interrupted  rota¬ 
tory  motion  of  a  permanent  (rigid)  dipole  is  represented  as  a  random 
walk  in  the  dipole-orientation  angle  space.  Collisions  are  considered 
to  be  instantaneous  events  (occurring  at  randomly  distributed  times) 
which,  while  preserving  the  dipole  orientation,  produce  maximal"*  random¬ 
ization  of  the  angular  velocity. 

Detailed  calculations  within  the  framework  of  this  model  have,  for 
the  most  part,  been  limited  to  the  domain  of  external  frequencies  small 
compared  to  the  natural  rotational  frequencies,  and  are  hence  not  rele¬ 
vant  to  the  present  problem.  The  extension  to  higher  frequencies  (of 

the  order  or  larger  than  the  rotational  frequencies)  was  carried  out  by 
6 

E.P.  Gross,  using  an  approach  based  on  a  Boltzmann  equation  in  orienta¬ 
tion-angle  space,  together  with  a  number  of  different  stochastic  models 
for  the  elementary  collision  process,  one  of  these  being  the  above 
described  Debye  model.  Instead  of  following  this  paper,  the  present 
authors  have  found  it  desirable  to  carry  out  an  equivalent  derivation, 
based  on  a  correlation-function  method This  derivation  is  presented 
in  Appendix  I.  Its  crucial  result  is  contained  in  equations  (1.24)  and 
(1.25),  which  show  that  for  external  frequencies  sufficiently  large 

g 

compared  to  thermal  rotational  frequencies  the  absorption  probability 
per  collision  (as  measured,  e.g.,  by  the  imaginary  part  of  polar¬ 
izability  is  essentially  identical  to  that  associated  with  the  collision- 
interrupted,  but  otherwise  free  motion  of  a  fictitious  charged  particle 
whose  charge,  e,  and  mass,  M,  are  related  to  the  actual  dipole  moment,  p, 
and  inertial  moment.  I  cf  the  rotator  by  the  formula 


2 

e 

M 


2 
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(1.1) 
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The  above  described  correspondence  is  in  fact  not  really  surprising. 

As  is  by  now  well  established  by  a  variety  of  treatments  on  pressure 

broadening,  the  absorption  intensity  at  a  given  frequency,  w,  is  deter- 

9 

mined  by  the  time-variations  in  the  electric  currents  on  a  time  scale 

'b  —  (or  less).  It  then  follows  immediately  that,  under  the  condition 
0) 

as  >>  tji  (postulated  above),  the  relevant  time  scale  is  so  small  that, 
rot 

within  it,  the  dipole,  so  to  speak,  "doesn't  know  that  it's  rotating. 

More  prosaically,  the  dipolar  current  (in  the  absence  of  collisions)  is 
essentially  constant  on  a  time  scale  ^l/tu;  the  significant,  and  meaningful 
current  variation  arises  solely  from  collisions. 

The  above  discussion  brings  us  naturally  to  the  consideration  of 
the  time-scale  of  the  collisions  themselves.  In  the  impact- formulation 
of  the  original  Debye  theory  (as  well  as  in  Gross'  work6  and  in  the 
appendix  of  the  present  paper),  this  time  scale  is  zero;  that  is,  a 
collision  event  causes  an  instantaneous  alteration  in  the  dipolar  cur¬ 
rent.  Such  an  assumption  is  appropriate  only  when  the  actual  time  of 
collision,  -  i.e. ,  the  time  required  for  the  angular  velocity  of  the 
dipole  to  change  —  is  small  compared  to  1/ to .  A  rough  estimate  of  this 
time  of  collision,  t  ,  is 


A 

•J-  %  - 

C  V 


(1.2) 


where  Ar  is  an  appropriately  defined  range  of  the  involved  forces,  and 
v  is  an  appropriate  velocity  of  traversal  of  the  force  field.  In  the 
at  ence  of  attractive  forces,  and  for  perturbers  whose  velocities 
(translational  or  rotational)  are  small  compared  to  the  rotational 
velocity  of  the  dipole  absorber, ^  a  reasonable  estimate  of  v  is 


v  n.  <2kT/Mred)- 


(1.3) 


where  M  ,  is  the  reduced  mass  for  rotational  motion;  for  e.g.,  HO 
red  L 

itself,  Mred*  »■„. 
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indicating  a  gross  violation  of  the  impact  condition. 

It  is  therefore  proposed  to  carry  out  a  calculation,  somewhat 
analogous  to  that  of  Chapter  I  of  Ref.  2,  in  which  one  computes  the 
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probability  that  a  charged  particle  incident  on  (and  reflected  from)  a 
potential  barrier  of  the  form 

V(x)  ■  V  e  aX  (1.6) 

o 

in  the  presence  of  an  electromagnetic  field  of  frequency,  u)  absorbs 

(o*  emits)  a  quantum  of  energy  fiu).  As  implied  by  the  form  of  (1.6)  the 

treatment  will  be  one -dimensional ,  with  x  representing  the  direction  of 

the  line  of  centers  of  the  colliding  atoms  (e.g.,  in  f^O  -  ^  collisions 

these  atoms  would  be  a  hydrogen  and  a  nitrogen  atom).  Of  course,  just 

as  in  Chapter  I  of  Ref.  2,  such  a  treatment  has  to  be  augmented  by  an 

appropriate-recipe  for  computing  effective  collision  diameters  (which  in 

a  one-dimensional  model  are  actually  infinite).  This  (rather  subsidiary) 

13 

problem  will  be  deferred  for  the  time  being. 
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2.  THE  MODEL  CALCULATION 

In  line  with  the  discussion  of  the  preceding  section,  let  us 
consider  the  following  problem:  to  solve  the  time-dependent  Schrodinger 
equation 

ifi  |-jf  =  Hip  (2.1) 

with  a  Hamiltonian  of  the  form 

H  =  H0  +  ^£-£  (2*2) 
ft2  2 

H  =  -  rr:  v  +  V(x)  •  (2.3) 

O  ZM 

In  these  expressions 

a  -  t  *'i"+ +  (2-4) 

is  the  vector  potential  associated  with  the  external  electromagnetic 
field  (it  is  assumed  to  be  "turned  on"  infinitely  slowly,  s  being  an 
infinitesimally  positive  number  ultimately  going  to  zero),  p  =  ^  grad^, 
and  Hq,  the  field-free  Hamiltonian,  is  given  as  a  sum  of  the  kinetic 
energy  of  the  fictitious  particle  (charge,  e,  and  mass,  M,  being 
determined  according  to  the  recipe  prescribed  by  the  Appendix  Eq.  1.27) 
and  the  potential  energy  V(x)  (assumed  to  be  one-dimensional  in  ac¬ 
cordance  with  the  discussion  of  Section  1).  In  the  absence  of  the 
external  electromagnetic  field,  (2.1)  has  solutions  of  the  form 

-iE . t/fi 

4»i(r,t)  =  ^(r)  e  ,  (2.5) 

where  4>^(£.)  satisfies  the  equation 

Ho4»i  =  Ei«p±  .  (2.6) 
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In  view  of  the  one-dimensional  form  of  the  potential,  <)K(r)  is  chosen  to 
have  the  form 


i[k  y+k  z] 

^(r)  =  <i1(x)  e  y  Z  (2.7) 

where  if)^(x)  obviously  satisfies  the  one-dimensional  equation 


1 

.  Lei 

[  2M 

£2  +  VW]  *  Eix  *i<x> 

(2.8) 

with 

Ei 

-  E,  + 
ix 

id  42  +  k2\  . 

2M  ly  z)  ’ 

(2.9) 

it  has  the  usual  asymptotic  form 


(x)  -»■  sin[kix  x  +  6^].  (2.10) 


Deferring  for  a  moment  the  complete  specification  of  the  boundary 
conditions  supplementary  to  (2.1),  let  us  seek  perturbed  solutions  of 
the  form 


-iE  t/ft 

<Hr,t)  =  ^(r)  e  +  ^ 


(2.11) 


where  if)  ^  is  proportional  to  the  external  field  (and  hence  arbitrarily 
small).  To  obtain  steady-state  solutions  (appropriate  to  the  physical 
situation  of  our  problem)  it  is  desirable  to  write  ip^  as 


♦j  ■  e 


|\(+)  „-ioit  A(»)  „+ii*it 


(+) 


where  <j>  satisfies  the  equation 


<f>v  /  e  +  <f>  7  e 

14 


(2.12) 


(E^-tfiu))  <p 


(+) 


H  <p 


P^Cr)  . 


(2.13) 
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Now  write  <p  ^  as 


e  &om  £  +1& 
iMw  -Rco  +  *£ 


(2.14) 


Substituting  (2.14)  into  (2.13),  one  finds  that  <j>  satisfies  the  equation 


(Ei+  1iw)*sc  -  Hq4>sc  =  ^  •  [Ho?  -  pHj^r) 


(2.15) 


e  &o- 

TujM  ftio  ‘fVP-pVH^r) 


•ox  3V(x)  , 


(2.16) 


(the  elimination  of  all  but  the  x-component  of  the  external  field  arising 
from  the  x-dependence  of  the  potential15). 

It  is  now  necessary  to  discuss  boundary  conditions.  As  will  be 
discussed  later  these  are  set  forth  most  expediently  for  <J>  [i.e.,  as 

supplements  to  Eq.  (2.16)]  rather  than  in  connection  with  the  starting 
Eq.  (2.1).  The  conditions  in  question  turn  out  to  be  ius.  that  <f> 

sc 

represent  an  outgoing  wave,  namely 


*  2^  e±[V  +  V1  A  e^ 

^SC  ®  SC 


(2.17) 


where  kfx,  the  "final"  x-component  of  wave-vector,  is  determined  by  the 
equation 


*2,2 
Ti  k. 
fx 


=  EfK  -  Ei+,i“-|i(ky  +  ^)  • 


(2.18) 


Here,  it  will  be  noticed  that  the  plane-wave  y,  z  dependence, 
characteristic  of  ^(r),  has  been  incorporated  explicitly  into  the 
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boundary  conditions.  In  fact,  at  this  point  it  is  highly  desirable  to 
factor  this  dependence  out  of  the  problem.  One  achieves  this  simply  by 
inserting 


_i[V 


k  z] 
z 


sc 


(2.19) 


into  (2.16).  One  obtains 


(x) 

E  ip  -  Hv  „ 

fx  sc  o  sc 


e^ox  3V  .  ,  n 

~2 —  Vx) 

w  m  3x  1 


(2.20) 


[H(x)  and  ih  (x)  being  defined  by  Eqs.  (2.8)  and  (2.7)]  whereas  the 
o  i 

boundary  condition,  (2.17)  reduces  to 


sc 


sc 


ik,  x 
fx 
e 


(2.21) 


With  the  y,z  dependence  now  completely  removed,  the  x-subscripts 
on  the  various  constants  will  henceforth  simply  be  dropped;  thus 


ix 


VEi 


ix 


krki 


(2.22a) 

(2.22b) 


(it  being  tacitly  understood  that  all  the  constants  now  refer  to  motion 
in  the  x-direction  alone). 

Before  going  on  to  solve  for  \p  (x) ,  some  brief  remarks  on  the 

SC  (+) 

physical  significance  of  the  decomposition  of  <£  as  given  by  the  r.h.s. 
of  (2.14),  are  in  order.  As  in  analogous  treatment  in  Ref.  2  (Chapter  1, 
bottom  of  page  11),  the  first  term  represents  the  amplitude  for  virtual 
absorption  of  a  quantum;  in  particular  one  notes,  upon  inserting  the 
asymptotic  form  of  iJj^(x)  [given  by  (2.7)]  into  (2.14),  that  the  asymptotic 
behavior  of  the  term  in  question  pertains  to  a  particle  having  a  kinetic 
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mmmmm 


■fi2k2 

energy  ,  characteristic  of  the  unperturbed  wave.  Indeed,  it 

would  even  exist  in  the  absence  of  any  scattering  potential;  in  fact, 
for  a  free  particle  its  coherent  superposition  with  the  unperturbed 
state  gives  just  the  classical  oscillatory  velocity  amplitude,  e  ^/iMco. 

By  way  of  contrast,  ^  [as  seen  explicitly  from  (2.21)]  describes 

s  c 

an  outgoing  particle  of  energy  =  E^+tioi.  This  fact  by  itself  is  suf¬ 
ficient  to  establish  it  as  the  amplitude  for  the  real  process  of  quantum 
absorption,  under  consideration  here;  it  will  therefore  alone  be  con¬ 
sidered  in  what  follows. 

From  these  remarks  it  follows  that  the  probability  of  an 

incident  particle  being  reflected  with  simultaneous  absorption  of  a 

quantum  of  electromagnetic  energy  is  proportional  to  A  2 ;  it  is  in 

16  SC 

fact  given  by  the  formula 


(2.23) 


For  the  calculation  of  Agc  one  proceeds  as  follows.  One  multiplies 
both  sides  of  (2.20)  by  that  solution,  i]»^(x),  of  the  homogeneous  equation 


H^f  =  Ef*f  (2.24) 

which  has  the  asymptotic  form 

<J>f(x)  sin(kfx  +  6f)  (2.25) 

and  integrates  from  minus  infinity^3  to  some  large  value  of  x  such  that 
(2.21)  and  (2.25)  are  valid.  One  then  obtains  the  successive  equalities 
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_ - £5.  1  ^(x ' ).  ^  (x’)ij;  (x')dx' 

<u2M  3x  '  Vx  'VX  )dX 


Ef  <Psc(x’)\pf(x')dx’ 


'X  fX 

^f(x,)Vsc(x,)dx'  "  KcVf  "  Wsc]dx' 
—  00  —00  ^  * 


‘  -  2M  |*,c(x)  to  +£  -  *f  A  *.c 


-  4  A  elkfX  to  sin(k£x  +  S£>  '  8ln<kfX  +  V  to  Asc  e 
2M  sc 

4  -i5f 
2M  Asc  6 


(2.26) 


the  next  to  the  last  equality  holding  by  virtue  of  footnote  17  and 
Gqs.  (2.21)  and  (2.25). 

Combining  (2.26)  and  (2.23),  one  has 


16ez£2 

^ox 


*f(5t)*iw£dx  2 


(2.27) 


where  the  limit  x  -*■  °°  has  been  taken  in  the  l.h.s.  of  (2.26),  the  limit 
being  well  defined  by  virtue  of  the  asymptotic  behavior  of  the  various 
factors  in  the  integrand. 

For  unpolarized  light,  the  standard  replacement 


g  2  1-2 

^ox  3  o 

is  tc  be  employed;  (2.27)  then  becomes 


16e2£  o  r  dv 
V)  -sror  J>(x)^i(x)^dx 


(2.28) 
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For  the  case  of  an  exponentially  repulsive  potential 


V  =  V  e 
o 


(2.29) 


the  integral  on  the  r.h.s.  of  (2.28)  may  be  evaluated  explicitly.  The 
procedure  is  very  similar  to  that  of  Ref.  2,  Chap.  I,  but  will  be  given 
here  briefly  in  the  interests  of  a  self-contained  presentation. 


The  solution  of 


E\p  = 


'  V  '  ['ll 


.  ,,  -ax  , 

TT  +  V„  e  'P 

dx*  o 


(2.30) 


(for  arbitrary  E)  with  the  asymptotic  behavior 


sin(kx  +6) 


(2.31) 


r  7  .<,1/2 

*<*)  -  [(2k/ an)  sin  -f-J  K2  ({) 


(2.32) 


where 


k  =  (2ME/fi2)1/2 


(2.33) 


-  (f) 


-ax/2 


(2.34) 


and  K  (?)  is  the  modified  Bessel  function  of  the  second  kind;  for  the 

2ik 

case  at  hand,  the  order  v  =  ■  ■  ■  is  imaginary. 


Inserting  (2.32)  and  (2.33)  into  (2.28)  one  obtains,  after  a  little 
algebra 


(+)  3tt 


£§hr  .In  sin  ^  |  1 2 

zMzoj4  a  a  1  1 


(2.35) 
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where 


X 


K2ik f/«  «>  K2ik. /a  «)5d5 
of  1 


(2.36) 


As  in  Chapter  I  of  Ref.  2,  it  is  our  good  fortune  to  find  the  above  integral 
1  P 

in  the  literature-  The  general  formula  (on  page  334  of  the  quoted 
reference)  is 

K^(ax)Kv(ax)xS  1dx  •=  2  r^(s+y+v)]  r[|-(s-y+v)] 

Jo 

x  r^(s+y-v)j  r^(s-y-v)]  ; 


a  >  o,(^e  s  >  |(5?ey|  +  |<^e  v|  .  (2.37) 


Here,  with  a  -  1,  s  =  2,  y  =  2ikf/ot,  v  =  21^/a,  the  conditions  for  the 
validity  of  (2.37)  are  obviously  fulfilled.  Using  also  the  well-known 
formula 


r(l+z)r(l-z) 


ZTT 


sinuz 

one  obtains,  without  further  ado 

1  (ir2/a2)  (k2-k2) 


k 


jinh^(kf+ki)j  sinh[a(kf"ki)] 


1 

2 


(it2/ a2)  (2Ma)/"ft) 


iinh^(kf+ki)J  sinh[^kf"ki)] 


7T  2 

7 


2Mqj/fi 


/2vkf\ 

cosh^ — - — -J  -  cosh 


(=9  ' 


u 


(2.38) 


(2.39) 
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Substituting  (2.39)  into  (2.35),  one  has 


P 


(+) 


167T2e2£2  sinh  (27rk.  /a)sinh (2Tik./a) 

222  Q  _  I 

3fi  a)  a  (Cosh  2?Tk^/a  -  cosh  2?rk^/a)^ 


(2.40) 


For  the  eventual  calculation  of  an  absorption  cross-section,  it 
is  desirable  to  develop  an  expression  for  the  net  energy,  AE  ,  which 
is  transferred  from  the  electromagnetic  field  to  the  charged  particle."^ 
To  obtain  AE^  it  would  appear  necessary  to  supplement  (2.40)  with  an 
explicit  expression  for  the  alternate  process  of  stimulated  emission 
[i.e.,  to  calculate  P._.].  This,  however,  can  be  avoided  by  appeal  to 
the  principle  of  microscopic  reversibility,  which  in  the  case  at  hand, 
states  that 


o 


P(+)(ki  ■“  kf)  =  P(-)(kf  ki)  * 


(2.41) 


One  now  introduces  this  relationship  into  the  general  expression  for 

AE  ,  which  reads 
0) 


AE  -  fiw 
0) 


‘■'o 


d^Kk^  P(+)(k1  -  kf)  - 


dkfI(kf)P(_}(kf  -  k±) 


(2.42) 


1} 
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where  e.g., 


Kk^dk.  = 


-E./kT 
kidki  e  1 

-E  /kT 
kidki  e  1 


dEi  -Ei/kT 

Tt"  6 


represents  the  probability  that  the  wave  vector 
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(2.43) 


and  where,  as  before, 


of  the  incident  particle  lie  between  k ^  and  k.-kik( , 

k.  and  k.  are  related  to  each  other  by  the  energy-conservation  require- 
f  1 

ment 

Ef  =  E1  ±  fico  *  (2.44) 


One  then  has  (upon  noting  also  from  (2.44)  that  k^dk^  -  k^dk^) 


AE  =  Tico 
0) 


0°  dE.  I  -E./kT  -E,/k.t\ 

‘h  1  -  l*(+)<Vkf> 


kT 


=  di 


4lVto/kT) 


»  -E./kT 


p(+)(ki>kf)dVkT 


(2.45) 


where,  for  (k£*kf) ,  one  inserts  the  r.h.s.  of  (2.40). 

The  stage  haa  now  been  reached  where,  in  order  to  make  further 
progress,  some  connection  is  needed  between  the  one-dimensional  treat¬ 
ment  of  this  section  and  the  actual  three-dimensional  situation.  The 
analogous  problem  arose  in  Chapter  I  of  Ref.  2  (pages  31  et  seq.),  and 
was  resolved  by  appeal  to  the  so-called  Taka', anagi  approximation  (for 
discussion  and  reference  to  Takayanagi's  paper,  consult  Ref.  2).  The 
principle  result  was  that  [as  intimated  above  in  footnote  (19)]  the 
•  ce  of  absorption  of  energy  per  polar  molecule,  W,  is  given  by  a 
relationship  of  the  form 


W 


AE  /t 
co  ef  f 


(2.46) 


where  1/ x  ,r,  the  "effective"  collision  rate,  is  to  be  determined 
ef  f 

primarily  by  classical  kinetic  conditions,  similar  to  (although  more 
complicated  than)  those  of  Ref.  2.  As  in  that  case,  however,  for  the 
purposes  of  obtaining  upper  estimates,  it  will  ultimately  be  assumed  that 


1/r 


eff 


>  1/t 


(2.47) 
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namely,  a  reliable  upper  bound  for  1/t^^  is  the  gas-kinetic  value,  1/t. 
Values  of  the  latter  quantity  are  readily  deduced  from  gas-kinetic  cross- 
section  data  given  in  the  literature. 

Combining  (2.46),  (2.45),  and  (2.40),  one  has 


W  = 


16ir2e2£^  g-lWkTj 


3a“ 


eff 


2 irk.  2irkf  -E  /kT 

sinh  -  sinh  — —  e  dE, /kT 

a  a  1' 


2irk, 


o  cosh 


-  cosh 


2irk.T 

ot 


(2.48) 


The  absorption  cross-section,  Q(w),  may  now  be  determined  simply 
by  dividing  (2.48)  by  the  incident  electromagnetic  flux, 

£  x  6o‘  0ne  has 


Q(w) 


.  -jftu/kT 
1-e _ 

hu 


1 

Teff 


2"k  2"kf  -E  /kT 

sinh  -jj-i  sinh  e  dEi/kT 

£cosh  2irk^/a-cosh  2iTki/aj 


(2.49) 


From  (2.49),  one  obtains  the  absorption  coefficient,  k (u>)  by  the 
standard  relation 


k(w)  =  N  0(b) 

a 

where  N  is  the  density  of  absorbing  molecules. 

3 


(2.50) 
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It  is  now  of  interest  to  discuss  two  special  cases,  (a)  that  for 
which  the  impact  approximation  is  valid  and  (b)  the  more  realistic 
situation  of  finite  duration  of  collision-time. 

A.  Impact  approximation ;  This  is  achieved  most  conveniently 
by  taking  the  limit  a  <*>  (i.e.,  zero  range  of  interaction 
potential.  In  this  case  one  replaces  the  hyperbolic 
functions  by  their  lowest-order  nonvanishing  terms  in 
their  Taylor  series,  obtaining  (after  some  algebra) 


Q(a)) 


4e2 

'  -WkTl 

1-e 

l  xl/2 

3irMco)2 

;  “ROJ  J 

■*0 

4e2 

Hio/kT 

e 

3ttMco)2 

2 

(x+f  ai/kT)1^2 

K  (1iu>/2kr) 


-x 

e  dx 


(2.51) 


c 


(the  integral  being  obtainable  in  terms  of  the  Modified 
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Bessel  Function  ).  In  the  classical  limit  (fiw  <<  kT)  (2.51) 
reduces  to 


4ez 

Qcl^  SttMco^ 


(2.52a) 


which,  apart  from  the  factor  4/3,  corresponds  to  the 

23 

classical  Drude  expression* 


In  the  "quantum"  limit,  'Tiu  >>  kT 

1/2 

4e2  j’ffkT  \ 

^qu  3irMca)^  l  2fiu)  I 


(2.52b) 


From  the  practical  standpoint,  the  principal  significance 
of  (2.52b)  is  that  its  order  of  magnitude  is  essentially 
the  same  as  that  of  the  classical  expression  (2.52a)  as 
long  as  fico/kT  is  not  too  large.  In  our  case,  with  -fiu)/kTV3, 
so  that  the  magnitudes  of  (2.52a)  and  (2.52b)  are  comparable. 
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B.  Finite  collision  times:  A  useful  simplification  of  (2.49) 
which  will  be  employed  here  is  the  replacement  of  the 
various  hyperbolic  functions  in  the  integrand  by  exponentials. 
For  the  justification  of  this  step,  let  us  digress  momen¬ 
tarily,  to  list  some  typical  magnitudes  of  l^,  kf,  and  a. 

First  of  all,  as  noted  in  Section  I  [text  subsequent  to 

8  -I 

Eq.  (1.51)],  a  representative  value  for  a  is  4x10  cm  . 

For  the  k's,  one  proceeds  from  the  general  expression 

k  =  ^  1/2  =  2.3  x  109(tt[E)1/2cm-1  (2.53) 

where  "Ml  is  mass  of  the  particle  in  units  of  hydrogen 
atom  mass  and  E  is  the  energy  expressed  in  e.v.  One  then 
notes  that,  even  for  the  extreme  case  of  iTl  =  1,  and  for 

g 

E  *=  kT  **.025,  one  has  k«»3.7  x  10  and,  e.g., 


2ifk  1  2irk  /  a 

sinh  - ‘t  e 

a  2 


-4nki/a  2vk^/a 

1-e  *  ~2  e 


1-e 


-11 


2irki  1 

so  that  the  replacement  sinh  ~ J  e 


involves  the  negligibly  small  error,  e 


-11 


Introducing,  then,  the  replacements  of  hyperbolic  functions 
by  appropriate  exponentials  in  (2.49),  one  has 


-E./kT 


8ire2 

f.  -Hco/kT  \ 
1-  e 

l 

l 

e 

dE  /kT 

3ca2 

[  11(0  I 

Teff 

4  sinh2 

0 

[z  <kf-ki)] 

2ire2 

L  -luo/kT ' 
1-e 

1 

-E  /kT 

e  1  dEi/kT 

3ca2 

t  Hw 

'  Tef f 

sinh2 

0 

2irai 

a(vf+  v±) 

(2.54) 


23 


where  v±  -  (2Ei/M)1/2  and  vf  «  (2/M)1/2 (E±+  ^)1/2  are  the 
velocities  incident  and  scattered  charged-particle, 
respectively . 

One  may  now  consider  two  subcases,  according  to  whether 
the  argument  of  the  hyperbolic  sine  in  (2.54)  is  small  or 
large  compared  to  unity. 

n\  - 2tuo —  «  1  In  this  case  the  hyperbolic  sine  is 

(vi+  V 

replaced  by  its  argument  and  one  has 


e2 

1 

bcirai^ 

\ 

Teff 

* 

e2 

1 

3'irca)'6 

\  tiu/kT  / 

Tef  f 

-E ,/kT 


(v^  vf)  e 


dE^kT 


x1/2  +  (x+tia)/kT)1/2 


-x 


(2.55) 


which  is  again  an  impact-type  expression  (somewhat 

24 

similar,  in  fact,  to  the  first  equality  of  (2.51). 


(2)  — - -  »  l.  In  this  case,  the  replacement 

v  a(v±+  vf) 

_  -  ,  1  (...) 
s  inh  (  ...  )  -*■  ~2  e 


is  appropriate,  and  one  has 


Q 


,00 


2 

•ne* 

,  -*WkTl 
1-e 

1 

3ca2 

|  -two 

Teff 

-E^kT 

e 


-2iro)/a(v  +vf) 

.  1  1  dE^kT 


dx 
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which,  with  this  introduction  of  a  new  integration 
variable,  y  =  E^/ftoi,  becomes 


Q 


,  -lWkT 

i-e 

kT 


e-bF(>’>  dy 

Jo 


where 


(2.56) 


and 


F(y) 

ay  +  1/2  1/2 
y  u+(y+l)lu 

(2.57a) 

_  ah_ 

a 

irkT  l  2M  1 

(2.57b) 

1/2 

b 

- 

a  rnii)  1 

(2.57c) 

The  integral  in  (2.56)  is  evaluated  in  Appendix  II;  the  result  is 
contained  in  Eqs.  (II. 7),  (II. 6),  and  (II. 8).  Combining  these  equations 
with  (2.56)  and  (2.57a,b,c),  one  has 


Q 


2 

•rre^ 

-fioj/kT' 

l-e 

1 

/4irkT\ 

<  *0)  J 

1  1  1  1 

3ca2 

kT 

T 

eff 

L».  V1  ”  ,»*+(. +iH 

o  S 

W  [yS+25's/2<V1)1/2D 


where  (cf.  Eq.  II. 4)  of  Appendix  II  and  (2.57b),  yg  is  given  implicitly 
by  the  relation 


ys/2(ys+1) 


1/2 


ys/2+(ys+1) 


1/2 


wkT 

[_M_V 

aft 

l  2tiaj  / 

1/2 


(2.59) 
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Values  for  a  and  v  =  u)/2irc,  for  various  values  of  y  „  are  given 
25  S 

in  Table  I . 


TABLE  I 


a 

vxlO  "^crn  1 

.1 

1.105 

3.75 

.15 

.826 

2.10 

.2 

.664 

1.36 

.25 

.553 

.943 

.3 

.476 

.695 

.4 

.367 

.415 

.5 

.299 

.276 

From  this  table  one  may  find  values  of  y  for  given  v  by  e.g. 

s 

graphical  interpolation.  It  would,  in  fact,  be  desirable  to  carry  out 
sufficient  numerical  calculations  so  as  to  obtain  a  plot  of  the  absorp¬ 
tion  cross-section,  Q,  vs.  external  wave-number,  v.  Due  to  limitations 
of  time,  this  elaboration  will  be  deferred  for  the  final  report.  How¬ 
ever,  in  the  next  section  (immediately  below),  a  numerical  upper-limit 
estimate  for  Q  (based  on  Eqs.  (2.58)  and  (2.59)  will  be  carried  out. 


3. 


DISCUSSION  OF  RESULTS  AND  POSSIBLE  EXTENSIONS 


The  first  order  of  business  in  this  section  is  the  comparison  of 

predictions  of  the  foregoing  theory  for  the  optical  absorption  coef- 

3 

ficient,  i.e.,  with  the  experimental  results  reported  by  Burch. 
Starting  from  the  relation 


k  =  NaQ  (3.1) 

where  N  ,  the  density  of  absorbing  molecules,  is  related  to  the  cor- 

responding  pressure,  p  (expressed  in  torr,  mm.  Hg)  via  the  formula 

a 

N  -  3.2  x  1016  p  cm"3  (3.2) 

a  a 

and  introducing  (2.58)  with  — - —  written  as 

Teff 


1 

Teff 


N  <  v> 

P 


^eff 


(3.3) 


(with  Np,  the  density  of  perturbing  molecules  given  by  the  analogue  of 
(3.2)  namely. 


N  =  3.2  x  10^  p  cm  3 

P  P 


(3.4) 


<v>  the  mean  relative  velocity  of  the  colliding  molecules,  and  Q 
and  "effective"  collision  cross-section,  one  has 


eff 


ire2N 


N 


eff  a  3a  ^kT 


<  v>  [~4kTl 

c  Liiii)  J 


1/2 


f"  exp^J . 


(3.5) 


where  from  (2.58) 


p"  = 


1-1/2 


Ju  +  -i_  +  _i _ l _ 

y  y  +1  1/2  ,  ...1/2 

s  ^s  yg  (ys+1) 
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\ 


mjjeaai 


and  the  argument  of  exp  /. . .  .^>  is 
given  by  (2.59),  or  Table  I. 


ftO) 

kT 


ys+2y 


1/2 

s 


(ys+1) 


y  being 
s 


Before  any  numerical  calculations  are  possible,  it  is  necessary 

to  determine  in  terms  of  knowa  physical  parameters.  Anticipating 

the  results  of  the  numerical  calculation  given  below,  let  us  bypass  this 

problem  by  adopting  the  philosophy  that  we  are  calculating  an  upper 
26 

bound  for  k.  In  line  with  this  approach,  we  write  down  an  upper- 

limit  estimate  for  Qg££  By  setting  it  equal  to  the  gas-kinetic  cross- 

section,  Q^.  For  the  latter  quantity,  a  rough  estimate  is  provided  by 

Laudolt-Bjornstein  (6th  Ed.)  Vol.  I,  p.  369,  where  a  value  for  the  gas- 

-8 

kinetic  diameter  o^,  of  H^O  —  namely  4.59  x  10  cm  —  is  quoted.  Intro¬ 
ducing  this  value  into  the  formula 


(appropriate  for  H, 0  -  H_0  collisions),  one  has 
i  z 


Qeff  ■  Qk  *  66  x  10  16cm2  . 


(3.6) 


8  —1 

Using  for  a  the  previously-cited  value  of  4  x  10  cm 
T  =  300°K  [together  with  formula  (3.2)],  one  finds 


and 


ire2N 


3a2kT 


2.  ■S 


1.2  x  10_6p  . 

P 


(3.7a) 


Taking  for  the  mean  relative  velocity  of  H„0  -  H_0  collisions  the 
27  L  1 

relation 


one  has 


<  v  > 


.76  x  10  cm/sec. 


<  v  > 
c 


2.5  x  10 


-6 


(3.7b) 


c 
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In  proceeding  further,  it  is  convenient  to  choose  for  the  external 

wave  number  the  value  v  =  934  cm  which  corresponds  (cf.  Table  I)  to 

3  -1 

y  =  .25  and  a  =  .553  (and  is  nevertheless  sufficiently  close  to  10  cm  ). 
s 

One  then  has 


(■*!?)  =  1,78  (3,7c) 

[F"]"1/2  =  .390  •  (3.7d) 


Inserting  these  results  [together  with  (3.4)  and  (3.6)  into  (3.5)], 
one  finds 


k  <  4.4  x  10  ^'pQp  exp{....}  .  (3.8) 

Turning  now  to  the  experimental  results,  Burch  concluded  that  in  an 
atmosphere-like  environment  with  a  l^O-vapor  pressure  of  15  torr  (mm  Hg) , 
self-broadening  (i.e,,  H20  -  H20  collisions)  constitutes  the  main  ab¬ 
sorption  agent.  In  xine  with  this  conclusion,  let  us  therefore  insert 

p  =  p  =15  torr 
a  *p 

into  (3.8),  thereby  obtaining 

k  *<"  10  exp{....}  .  (3.9) 

Deferring  for  the  moment  the  numerical  evaluation  of  exponential 
factor,  we  note  that,  under  the  cited  conditions,  Burch  obtains  .88 
transmission  of  10  micron  light  through  a  path  of  1  km;  this  cor¬ 
responds  to  an  absorption  coefficient 

k  =  -10"5  log  .88  ‘S’  1.3  x  10"6cm_1  (3.10a) 

ex 
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which  is  one  order  of  magnitude  larger  than  the  pre— exponential  factor 
in  (3.9). 

Turning  finally  to  the  exponential  factor,  one  finds  [from  the 
expression  for  its  argument  given  in  the  text  immediately  subsequent  to 
(3.5)]  that  (with  WkT  =  2Trficv/kT  =  3.98) 

exp{....}  -  exp{-(3.98) (1.37)}  =  exp{-5,A5}  =  .A3xl0-2  .  (3.10b) 

It  is  thus  seen  that  the  theoretical  value  of  k  is  smaller  than  Kexp 
by  at  least  three  orders  of  magnitude! 

In  seeking  a  road  out  of  this  impasse,  one  should  note  that  the 
primary  "culprit"  is  the  exponential  factor,  accounting,  as  it  does, 
for  two  out  of  the  three  orders  of  magnitude  of  discrepancy.  Physically, 
this  factor  represents  the  effect  of  finiteness  of  collision  time 
in  particular  the  fact  that  Tc  is  so  large  that  the  condition  for  the 
validity  of  impact- type  theory  (namely,  u>Tc  <  1)  is  grossly  violated. 

In  this  connection,  it  may  be  helpful  to  recall  the  remarks  of 
the  second  paragraph  of  Chapter  I  of  Ref.  2  (albeit  somewhat  rephrased 
for  application  to  the  present  problem).  Three  alternative  physical 
regimes  were  therein  described. 
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1.  The  domain  of  external  frequencies  obeying  the  condition 

air  <  Is  This  is  the  regime  of  the  impact  theory, 
c 

2.  For  frequencies  such  that  uit^  >  1,  the  spectral  profilo 
will  follow  the  predictions  of  the  so-called  statistical 
theory.  In  this  theory,  one  focuses  attention  on  the 
collision-induced  perturbation  of  the  instantaneous 
transition  frequency.  Two  possibilities  arise.  If 

a  collision  gives  rise  to  perturbation  of  the  system 
(during  the  transient  existence  of  the  collision  complex 
H20  -  H20),  such  that  a  perturbed  transition  frequency 


c 
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momentarily  coincides  with  the  external  frequency,  u>,  the 
spectral  intensity,  I  (to),  will  be  proportional  to  the 
occurrence- probability  of  such  a  perturbation  —  i.e., 
the  probability  that  any  one  of  the  perturbed  (rotational) 
transition  frequencies  lie  within  a  unit  different ial 
frequency  range  about  the  spectral  frequency,  to. 

3.  If,  on  the  other  hand,  there  are  no  collisions  which 
give  rise  to  perturbed  transition  frequencies  in  the 
vicinity  of  the  spectral  frequercy,  to,  one  may  expect 
the  spectral  intensity,  I  (to),  to  decrease  exponentially 
with  increasing  frequency-discrepancy  (said  frequency- 
discrepancy  being  the  difference  between  to  and  the  closest 
perturbed  (rotational)  transition  frequency). 

With  respect  to  the  theory  of  Section  II  of  this  paper,  we  have 
already  seen  that,  under  the  given  physical  situation  of  interest  to 
us,  the  condition  totc  <  1  does  not  apply,  so  that  the  impact-type 
frequency  distribution  (case  1)  does  not  occur. 

As  far  as  the  second  of  the  above-listed  three  cases  is  concerned  — 

the  theory  contains  no  effective  mechanism  for  rotational  perturbations 

of  the  magnitude  required  to  encompass  frequencies  in  the  domain  of 
3  -1 

interest  (10  cm  ).  The  frequency  variation  of  the  absorption  coef¬ 
ficient  therefore  falls  into  the  third  category  —  as  is  in  fact 
evidenced  by  the  appearance  of  the  exponential  factor  on  the  r.h.s.  of 
Eq.  (3.5)  (the  negative  argument  of  which  is  readily  shown  from  Eq.  (2.58) 
and  Table  I  to  be  a  numerically  increasing  function  of  frequency). 

At  this  point  a  qualifying  remark  is  in  order.  In  a  strict  mathe¬ 
matical  sense  the  change  in  rotational  motion  which  arises  via  inter¬ 
action  of  the  rotating  dipole  with  the  (assumed  exponentially  repulsive) 
potential  of  the  colliding  particle  —  taking  place  as  it  does  in  a 
finite  time,  ic  —  is  characterized  by  a  frequency-occurrence  distribution 


31 


3-1 

which  certainly  includes  the  spectral  frequency  of  interest  (v  =  10  cm  ) 
(otherwise  the  absorption  coefficient  at  that  frequency  would  be  zero). 

In  fact,  from  this  strictly  mathematical  point  of  view,  there  is  no 
distinction  between  cases  2  and  3.  However,  especially  under  the  con¬ 
ditions  oil  >>  1,  it  is,  in  the  opinion  of  the  present  author,  physically 
c 

meaningful  to  ask  whether,  within  some  sub-interval  of  the  collision 

time,  x  ,  the  rotational  motion  may  be  characterized,  at  least  momen- 

c  3-1 

tarily,  by  a  frequency  in  the  vicinity  of  10  cm  .  Such  a  possibility 

will  in  fact  be  suggested  at  the  end  of  this  section;  although  detailed 

calculations  have  not  been  performed,  it  is  (hopefully)  anticipated  that 

the  suggestion  will  provide  the  key  to  the  understanding  of  Burch's 

experimental  findings. 

However,  to  proceed  systematically,  let  us  first  investigate  one 
other  possibility  of  increasing  the  absorption  cross-section  —  a 
possibility  suggested  by  comparison  of  (2.58)  with  the  previously 
obtained  impact  expression  (2.55).  Such  comparison  —  taken  in  con¬ 
junction  with  the  text  immediately  subsequent  to  Eq.  (3.10)  — 
immediately  poses  the  question  as  to  whether  some  mechanism  can  be 
found  to  diminish  the  argument  of  the  hyperbolic  sine  in  Eq.  (2.54) 
to  the  point  where  sub-case  1  of  the  previous  section  obtains,  i.e., 
where 


2  -TO) 

a(vi+vf) 


<  1  . 


(3.11) 


Such  a  mechanism  is,  at  least  in  principal,  provided  by  the  existence  of 
attractive  interactions;  the  question  as  to  whether  it  is  actually  ade¬ 
quate  will  now  be  studied. 


In  order  to  simplify  the  discussion,  let  us  introduce  the 

approximation  of  considering  the  attractive  interaction,  V^(x)  (which, 

in  practice  would  vary  as  some  inverse  power  of  x.. )  to  be  constant, -V  , 
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—  equal  in  magnitude  to  its  value  at  the  classical  turning  point. 
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Just  as  in  the  somewhat  analogous  Takayanagi  approximation,  discussed  in 

<>  Ref.  (2),  (Chapter  I,  p.  31-32),  this  approximation  should  be  valid  —  at 

least  for  semiquantitative  estimates  —  when  V  (x)  varies  sufficiently 

& 

slowly  so  that 


ki  I  dbc  log  Va(x)  I  >  1  (3.12a) 

kf  |  log  Vfl(x)|  >  1  .  (3.12b) 

Since  this  condition  is  already  on  the  borderline  of  being  obeyed  by  the 
relatively  rapidly  varying  repulsive  potential,  it  will  here  be  assumed 
to  be  applicable  without  further  ado. 

With  the  assumption  of  constant  attractive  potential  the  analysis 
becomes  rather  simple,  in  that  the  velocities  v^  and  v^  in  (2.54)  are  now 
given  by  the  formulae 


v±  =  (2/M)1/2  (E±  +  Va)~/2 
vf  =  (2/M) 1/2  U  r  V  +  bui)1/2  , 

X  m  3. 
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with  retaining  its  former  significance  as  the  initial  energy  at 

x  •*  00 .  _ 


(3.13a) 

(3.13b) 


1 


Introducing  (3.13a)  and  (3.13b)  into  Eq.  (2.54)  —  together  with  the 
additional  assumption  (to  be  verified  below)  that  for  the  case  of  interest 
Va  is  sufficiently  larger  than  to  permit  us  to  neglect  the  latter 
quantity  in  (3.13a)  and  (3.13b),  the  integral  in  (2.54)  is  immediately 
evalued,  and  one  obtains 


(3.14) 


One  sees  that,  provided  that  V  is  sufficiently  large  i.e.,  large 

Sl 

enough  so  that  the  condition 


2TTa)/a 


(2/M) 1/2  [vy2+(Va+ttw)1/2] 


<  1 


(3.15) 


32 

is  obeyed,  impact- type  behavior  will  result,  with  Q  being  given  by 

,2 

e  1-e - -  - 2 - i_  .  (3.16) 


2  ,  -TWkT  M'W'2**)] 

i  1-e  l  a _ a _ J 


3ttMo)2c 


eff 


If,  on  the  other  hand,  (3.15)  is  not  fulfilled,  (3.14)  becomes 


%  2 ire 


3ca 


1-e 


-WkT 


■hoi 


eff 


exp  (  - 


4iTco/a 


(2/M)1/2[vy2+(Va+ha>)1/2  ] 


(3.17) 


In  order  to  proceed  further  it  is  now  necessary  to  put  numbers  into 
(3.15).  Inserting  values  for  the  physical  quantities  which  have  already 
been  employed,  viz: 

/  ,n8  -1 

a  *  4  x  10  cm 

v  -  934  cm  ^ 

and  M  =  2 

into  (3.15),  one  obtains 

V  +  (V  +  ■hw)^2  >  2.9  ev  .  (3.18) 

a  a 


o 
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The  required  value  of  Vfl  is  then  seen  to  be  *  (2.9)2-  8.4  e.v.  Such  a 
magnitude  is  impossibly  large. 

A  reasonable  estimate  of  the  magnitude  of  V  to  be  expected  in 
collisions  may  be  obtained  by  using  force-constant  data  derived 
from  experimental  values  of  the  virial  coefficient  in  H^O  vapor.  A  fairly 
exhaustive  discussion  of  various  representations  of  the  interaction- 
potential  between  two  water  molecules  -  with  associated  numerical  values  — 
is  contained  in  the  book,  Molecular  Theory  of  Gases  and  Liquids  by 
Hirschfelder ,  Curtiss,  and  Bird  (HCB) .  For  our  purposes,  one  of  the 
simpler  forms,  the  so-called  Stockmayer  potential  [HCB,  page  211,  Eq. 
(1.3-33)], 


Hi‘H2 


3(Mi*r)(y2*r) 


J 


j^2  cos  0X  cos  02  +  sin  ei  sin  02 


cos 


(3.19) 


will  be  employed.  In  this  form,  the  interaction  between  two  identical  polar 
molecules  is  expressed  as  a  sum  of  a  Lennard- Jones  potential  (r-intermolecular 
distance)  and  a  dipole-dipole  term,  in  which  the  charge  distributions  of 
the  individual  molecules  are  represented  as  point  dipoles  (0^  <j>,  and  e2, 

<t>2  being  the  polar  coordinates  of  the  dipoles  with  polar  axis  parallel  to 
the  intermolecular  vector,  r) .  The  values  of  the  constants  e,  and  a  for 
the  H20-H20  interaction  are  listed  in  Table  (3.10-1)  (page  216  of  HCB). 

They  are 


e/k  -  380°K 


(3.20a) 


a  -  2.65A 


(3.20b) 
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which,  together  with  the  known  dipole  moment 


-18 

^  q  =  1.84  x  10  e.s.u.  (3.20c) 


provides  us  with  a  concrete  expression  for  V(r). 

The  procedure  to  be  used  here  for  obtaining  an  upper  estimate  for 
V  from  (3.19)  is  as  follows.  One  maximizes  the  attractive  interactions 

8l 

by  taking  the  dipole  moments  and  ^  both  parallel  to  the  intermolecular 
separations  (i.e.,  0^  =  62  *  0) ,  (3.19)  thereby  reducing  to 


V  = 


(3.20) 


2  3 

One  then  identifies  V  with  the  values  of  oty  /r  ,  where  r  is  the  closest 

a  '  o 

distance  of  approach  of  the  two  molecules.  An  approximate  condition  for 

33 

the  determination  of  r  is 

o 


V  (rQ)  -  0 


(3.21) 


which  with  the  introduction  of  a  new  variable 

x  =  (o/r) 
o  o 


gives 


0  -  V(x  ) 
o 


(3.22) 


(3.23) 


Discarding  the  trivial  solution,  xq  =  «,  one  has  then  to  solve  the 
cubic  equation 
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o 


(3.24) 


From  (3.20a),  (3.20b),  (3.20c),  one  determines  the  numerical  value  of 
the  constant  term  in  (3.24)  to  be 


3 

2c  e 


A  sufficiently  accurate  solution  of  (3.24)  and  (3.25)  is 


(3.25) 


x  =  1.48 
o 


(3.26) 


which,  when  substituted  into  the  last  term  of  (3.23),  gives 


v(r,^x)  =  2_l  _  x  =o  5.4  x  10  erg  =  .34  ev  . 


3  3  o 

r  a 


(3.27) 


Comparing  (3.27)  with  the  value  of  V  required  to  satisfy  (3.18) 

i\,  ® 

—  namely  V&  >  8.4  ev  —  one  sees  that  (3.17),  rather  than  the  impact  ex¬ 
pression,  (3.16),  is  the  appropriate  formula  for  Q.  Let  us  therefore 
combine  (3.16)  with  (3.1)  and  (3.3)  to  obtain  the  following  expression  for 
the  absorption  coefficient 


K 


Qeff 


N 


a 


8tre2 

,  -hwkT 

1-e 

<v>  f 

4iTa)/a  **1 

3a2 

i 

‘  i'  (2/M)1'2 

V1/2  +  (V  +tla>)1/2 

3  Si  m 

J 

(3.28) 


Upon  inserting  the  numbers  given  by  (3.21),  (3.7a),  (3.6),  (3.7b)  and 


*  V 

%  i 
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(3.27)  into  (3.28),  and  utilizing  the  (appropriate)  formula  M  =  214^, 

as  well  as  a  =  4  x  10^  cm  \  v  =  943  cm  \  and  T  =  300  K,  one  has  (with 

,  .  ,  -haj/kT 

the  neglect  of  e  , 


K  ^ 


1.06  x  10  ^  PaPp  exp  ^-4.6^ 
2.5  x  10  ^  exp  ^-4.6^  cm  ^ 


-1 

cm 


(3.29a) 


(the  last  number  being  obtained  from  p  ■  p  "15  torr,  in  conformance 

a  P  -4.6 

with  Burch's  experimental  conditions).  Finally,  with  i  *  =  .01,  one 

has 


k  =  2.5  x  10  ^  cm  ^  (3.29b) 

which  exceeds  the  previous  upper  limit  estimate,  given  by  (3.9)  and 
(3.10),  by  a  factor  of  2. 5/. 43  =  5.8,  but  still  falls  short  of  the  ex¬ 
perimental  number,  given  by  (3.10),  by  a  factor 

k  / k  .  >  520  (3.30) 

ex  th 

i.e.,  still  close  to  three  orders  of  magnitude'.  It  thus  appears  that  a 

substantially  new  physical  ingredient  has  to  be  injected  into  the  theory. 

In  what  follows,  one  such  possibility,  which  appears  to  be  hopeful,  will 

34 

be  discussed  briefly. 

The  new  physical  ingredient  is  in  fact  contained  in  Ec}.  (3.19) 

—  namely,  in  its  dependence  on  the  angular  orientations  of  the  dipoles 
(as  given  by  the  second  term).  Preliminary  estimates  indicate  that,  at 
intermolecular  distances  of  the  order  of  rQ  "  c/xq  %  2.2A  (as  given 
by  (3.26)  and  (3.20b)),  the  rotational  motion  is  sufficiently  hindered 
so  as  to  be  converted  into  libration  about  the  potential  minimum 
(0^  3=  a  o).  According  to  (rather  crude)  preliminary  estimates,  the 


C 
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associated  librational  frequencies  are  in  a  domain  v  500  cm  .  In  view 

4#’  of  the  fact  that,  at  the  distance  of  closest  approach  the  dipole-dipole 

approximation  is  certainly  inaccurate  for  quantitative  estimates  of 

35 

librational  frequencies,  this  estimate  should,  in  the  opinion  of  the 
present  authors,  be  regarded  as  promising* 

In  this  connection,  some  remarks  concerning  the  absorption  spectrum 

3  6 

of  H^O  in  its  condensed  phases  are  in  order.  It  is  known  that  whereas 
the  vibrational  spectrum  of  1^0  vapor  is  largely  retained  in  both  water 
and  ice  (with  relatively  small  frequency  shifts) ,  the  rotational  spectrum 
is  wiped  out*  In  its  place,  there  appears  a  broad  absorption  band,  whose 
maximum  varies  from  700  cm  ^  (liquid)  to  850  cm  ^  (ice  at  -9°C) .  This 
band, usually  designated  as  a  "librational"  band,  is  considered  from 
hindered  rotational  motion  —  the  latter  being  associated  with  H-bond 
formation*  The  H-bond  in  turn  is  generally  considered  to  owe  its 
existence  to  the  primarily  electrostatic  attraction  between  one  of  the 
(positively-charged)  H-atoms  of  a  given  ^0  molecule  and  the  (negatively- 
charged)  0-atom  of  a  neighboring  ^0  molecule* 


Apart  from  quantitative  detail,  this  picture  is  essentially 

identical  to  that  suggested  in  this  section*  Namely,  even  in  the  gas 

phase,  at  sufficiently  close  distances  of  approach  of  two  colliding  ^0 

molecules,  a  temporary  dimer  (^O^  is  formed;  it  is  eminently  reasonable 

that,  during  its  (admittedly  momentary)  existence,  this  dimer  possesses 

librational  (hindered-rotational)  frequencies  in  the  spectral  range  of 
3  -1 

interest  (10  cm  )  <,  The  principal  problem  is  then  to  compute  the 

occurrence-probability  of  a  unit  frequency  range  about  an  arbitrary  spectral 

37 

frequency;  according  to  the  statistical  theory  of  collision  broadening 
a  knowledge  of  this  o  currence-probability  will  lead  directly  to  the 
absorption  coefficient  for  said  spectral  frequency* 


*  V 
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APPENDIX  I 

In  this  appendix,  the  relationship  between  the  real  problem  of 
absorption  associated  with  the  collision-interrupted  motion  of  a  rotating 
dipole  and  the  text-model  of  absorption  by  a  charged  particle  undergoing 
collisions  with  fixed  scattering  centers  will  be  developed. 

The  basic  starting  point  of  the  treatment  is  the  Kramers-Heisenberg 
formula  for  the  polarizability,  a,  of  a  (rigid)  rotating  dipole,  inter¬ 
acting  with  a  fixed  (random)  array  of  scattering  centers  (generalized  to 
include  the  dissipative  and  finite  temperature  effects) .  This  expression 
reads 


a 


nm 


<  n|bz|m> 


2 


1 

E  -E  -rtlu+ihs 
n  m 


+  E  -E  -ftco-itis 
n  m  J 

Here  3  -  1/kT  (with  k  and  T  respectively  Boltzmann's  constant  and 

absolute  temperature) ,  Z  is  the  partition  function  of  the  system  En  is 

any  one  of  the  (exact)  eigen-energies  of  the  total  system,<n|y  |m> 
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is  a  typical  matrix  element  of  e.g,,  the  z-component  of  the  dipole 

moment,  w  is  the  external  frequency  and  s  an  infinitesimal  positive 
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number. 

Of  specific  interest  here  is  the  imaginary  component;  of  polari¬ 
zability,  a^.  This  quantity  is  gotten  by  taking  the  imaginary  part  of 
the  r.h.s.  of  (1.1);  with  the  aid  of  the  standard  relationship 
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it  takes  the  form 
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(1.2) 


wherein  the  device  of  oummation-index  interchange  haa  been  utilized  re 
peatedly,  and  where  cognizance  has  been  taken  of  the  relationship 

E  -  E  +  hu 
m  n 


in  writing  the  third  equality. 

One  now  expediently  introduces  the  Dirichlet  representation 

s<*>  -27  /  •**  it 

—00 


-*1 1 

*  > 


for  the  delta  function  into  (1.2);  the  standard  rules  of  matrix  multi¬ 
plication  may  then  be  utilised  to  rewrite  (1.2)  as 


ctj  -  —  tanh 


6hw  f 

2  J 


cos  wt<p  (t)  y  (o)>  dt 
z  z 


(1.3) 
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where 


1  _  f  -SH  iHt/fi  -iHt/fi 

2  Tr  <  e  e  Mze  y 


(1.4) 


is  the  dipole-dipole  "response"  function,  H  being  the  Hamiltonian 
operator  of  the  total  system. 

In  this  appendix,  the  discussion  will  be  restricted  to  classical 
conditions  —  speef  ically^ 


fttD  <  kT 


(I. 5a) 


lito 


rotational 


<  kT 


(I. 5b) 


(1.3)  may  then  be  approximated  by 

cos  wt<  y  (t)y  (o)>dt  .  (1.6) 

z  z 

—00 


In  order  to  proceed  further,  specific  assumptions  concerning  the 
dynamics  of  the  system  have  to  be  introduced.  In  this  treatment,  said 
dynamics  will  be  described  by  a  stochastic  model,  of  the  general  type 
developed  in  the  theory  of  Brownian  motion.  Specifically,  it  will  be 
assumed  that 

(a)  in  the  absence  of  collisions,  the  dipole  undergoes  simple 
rotatory  motion  (characteristic  e.g.,  of  a  simple  rigid- 
rotator  possessing  no  internal  angular  momentum)  and 

(b)  at  certain  "impact"  times,  t  ,  the  vector  angular  velocity, 

changes  abruptly  to  a  new  value,  however,  the 


4i 

orientation  of  the  dipole,  y^/y,  remains  unchanged.  It 
is  further  assumed  that  the  direction  of  apart  from 

being  perpendicular  to  the  rotator-axis  (and  hence  to  jj^)  is 
otherwise  randomly  arbitrary. 


On  the  basis  of  this  model,  the  calculation  of  the  response  function 
<%z(t)  yz(o)^>  may  be  carried  out  as  follows.  In  the  time  interval  be¬ 
tween  zero  and  t  (here  considered  positive  for  the  sake  of  definiteness; 
later,  it  will  be  established  that  the  response  function  is  an  even 
function  of  t)  there  will  in  general  occur  a  certain  number  of  impacts  at 
times  tp  t2»  ....  tn_^;  *n  each  of  the  n  intervals  in-between  collisions, 

the  motion  is  that  of  pure  rotation.  Denoting  by  jj,  the  value  of  n  at 
til 

the  i  collision  time,  t^,  one  has 


U.(t)  -  cos  oj1  (t-ti) 


gj+1  X  Hj 

“i+l 


sin  (t-t.) 


(  ) 


Using  the  above  stated  assumption  concerning  randomness  of 

direction  of  the  angu1 ar-velocity  vector,  m.  . ,  in  a  plane  perpendicular 
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to  one  has  immediately  (upon  replacing  i  by  i-1  for  convenience) 


<u(t')>  -<j!1.1><co8  w1(t,-t)>;  (t±1  S  t  5  t±;  tQ  -  0).  (1.7) 


In  particular 

<£1>=<Hi_1><cos  mi  (t1-t1_1>  =<yi_1  cos  U±  0±> 

where  0^  ■  t^  -  t^  ^  the  time  interval  between  the  i'th  and  the  (i-l)^ 
collision;  the  various  6^  are  assumed  to  be  distributed  according  to  the 
standard  relationship 
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w(ei)dei  =  e 


de^t  . 


<I-&) 


-6 


For  the  final  interval  between  the  times  of  last  impact,  tn_^ 
time  t,  one  has 


<£(t)>  =<VlXC°S  Wn6n> 


with  0  *  t  -  t  n ,  the  occurrence  probability  of  0  (i.e.,  that  there 

n  n-1 

be  no  collision  in  time  6  )  being  given  by 


W 


n 


(I. 10) 


Putting  all  these  results  together  one  has  for  the  contribution  to 
y  (t)  y  (0)  due  to  all  processes  involving  n  collisions  between 

time  zero  and  t. 


<u,(t)  V2<0)> 
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o  o 


n 

n 

i=l 


cos  U) 


(1.11) 


On  the  r.h.s.  of  this  expression,  the  factor,  3,  in  the  denominator 
arises  from  taking  an  average  overall  initial  orientations  (at  t=0) ,  and 
the  delta  function  guarantees  that  the  sum  of  all  the  intervals,  0^,  is 
equal  to  t.  Finally,  the  superscript  "(+)"  on  the  l.h.s.  depicts  the 
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fact  that  the  relationship  is  valid  only  for  positive  time;  otherwise  the 
r.h.s.  is  manifestly  zero. 

With  respect  to  negative  times,  it  may  be  noted  that,  since  the 
correlation  function  depends  only  on  relative  time,  i.e., 

one  has 

<^yz(-t)uz(0)^>  =  <^yz(0)yz(t)^>  .  (I. 12a) 

Now,  in  a  classical  theory,  pz(t)  and  yz(0)  are  c-numbers,  so  that  the 
order  in  which  they  are  written  is  immaterial.  It  then  follows  immediately 
that,  for  t  >  0 

<yz(_t)yz(0)^>  “^yz(t*yz(0^ 
or 

<yz(t)yZ(°)>  =  ^'Jz(lt^yz(0^+)  (1.13) 

the  r.h.s.  of  (1.13)  being  given  by  (I. 11). 

With  the  aid  of  (1.13),  one  may  rewrite  the  basic  relationship 
(1.6)  for  as 


a.  =  gw  /  cos  wt  <u  (t)u  (0)>  dt  . 


(1.14) 


Since,  now,  only  positive  times  are  involved,  one  may  forthwith  introduce 


n 
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MH 


(I. 11)  into  (1.14);  a  subsequent  time  integration  then  yields  (upon  also 
summing  over  all  interval  numbers,  n,  from  one  to  infinity). 
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(1.15) 
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The  braces  signifying  the  remaining  average  over  each  of  the  w.. 

Concerning  this  average,  it  is  undoubtedly  most  appropriate  to 
assume  (as  in  fact  done  in  Gross's  paper^)  that  the  uk  are  each  thermally 
distributed.  However,  in  the  interest  of  transparency,  the  simpler  pro¬ 
cedure  of  taking  all  the.  uk  equal  to  one  and  the  same  angular  speed,  u>r 
given  eventually  by  an  appropriate  thermal  average  will  be  adopted. 


Introducing,  then,  the  replacement  ■*  u>r,  and  carrying  out  the 
sum  over  n,  one  obtains,  after  a  little  algebra 
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(1.16), 


or,  more  explicitly, 
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(1.17) 
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For  the  case  in  which 
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(1.18) 
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(which  applies  to  our  problem),  and  for  to  not  too  close  to  id  (i.e.,  in 
the  far  wings  of  rotational  lines),  (1.17)  reduces  to 
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(1.19) 


At  this  point,  it  is  desirable  to  take  the  thermal  average  over  id  . 
For  the  case  of  a  simple  (linear)  rotator,  the  weight  function  is 

-8I(d2/2 

to  e 


so  that 


/  2\_  _2_ 
^“r^  61 

Or>  =  2<(d2> 


(I. 20a) 


(I.20b) 


and  (1,19)  becomes 
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(1.21) 


^  somewhat  more  useful  form  of  (1.21)  is  achieved  by  expressing 
in  terms  of  the  characteristic  rotational  constant, 


From  (1.22)  and  (I.20a),  one  has 


(1.22) 
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(1.23) 


where  u>  has  been  expressed  in  wave  numbers 


v  =  u/2ttc  . 


Inserting  (1.23)  into  (1.21),  one  has 
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It  is  now  of  irterest  to  estimate  the  magnitude  of  the  second  term  in 
the  square  bracket  for  the  case  of  interest  to  us.  Here  the  relevant 
numbers  are  v  =  1000  cm  ^ ,  kT/2Trttc  'v  300  cm  ^  and^  B  30  cm  \  Insert¬ 
ing  them  into  (1.24)  gives  for  the  second  term  e  magnitude  %  0.15.  Let 
us  now  note  that  the  whole  effect  of  rotational  motion  is  contained  in 
this  tern).  The  above  numerology  exhibits  that  this  motion,  although  non- 
negligible  is  of  secondary  importance  in  determining  the  absorption  at 
10  microns.  In  particular  it  justifies  the  text  procedure  of  replacing 
rotatory  motion  of  the  dipole  by  linear  motion  of  a  fictitious  charged 
particle. 

A  final  point  concerns  the  values  of  charge  and  mass  to  be  used 
for  this  fictitious  particle.  In  the  case  of  a  diatomic  rotator,  one  may 
express  y  and  I  in  terms  of  internuclear  separation,  R,  as 


(1.21)  then  becomes 
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(1.25) 


The  corresponding  expression  for  the  real  part  0  of  the  complex 

R 

"conductance",  a  -  -ioxx  (i.e.,  o _.=  oi  .)  reads 
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(1.26) 


which,  apart  from  the  correction  term  in  the  square  bracket  (discussed 
above)  and  the  factor  2/3  (reflecting  the  fact  that  the  number  of  degrees 
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of  rotational  motion  is  two,  as  compared  with  three  for  translation)  is 
just  the  well-known  Drude  formula  for  absorption  associated  with  the 
(impact)  collision- interrupted  motion  of  a  free  particle  of  ch?rge  e 
and  mass  M  ed-  It  is  thus  clear  that  whenever  the  ratio  e  /Mred  occurs 
in  text  equations44,  the  appropriate  replacement  recipe  is: 


(I 


.27) 
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APPENDIX  II 


A  f 


In  this  appendix,  the  integral  in  the  text  Eq.  (2.56)  will  be 
evaluated  by  the  method  of  steepest  descents.  The  integral  in  question 
is  of  the  form 


& 


dy 


(II. 1) 


o 


where 


F(y)  =  ay  +  l/[y1/2  +  (y+l)1/2]  (II. 2) 

(the  dimensionless  parameters,  a  and  b,  being  defined  by  the  text  Eqs. 
(2.57b)  and  (2.57c),  respectively). 

The  first  step  is  to  locate  the  saddle-points,  as  defined  by  the 
equation 

F'(y)=0.  (II. 3) 


Differentiating  (II. 2),  one  has 


F'(y) 


v-1/2  +  (v+l)-1/2 
'  2[yW*  +  (y+l)1/2]2 


which,  after  a  few  algebraic  manipulations,  takes  the  form 


F'(y) 


(II. 4) 
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Introducing  (II. 4)  into  (II. 3),  one  obtains  the  saddle-point  equation 


2yJ/2  <ys+D1/2 


1/2  .  t  N 1/2 

ys  +  (ys+D 


(II. 5) 


(the  subscript,  s,  denoting  the  saddle  point  values4^  of  y) . 

The  usual  procedure  is  to  solve  for  y  in  terms  of  the  (physically 

s 

variable)  parameter,  a.  In  view  of  the  complicated  algebraic  form  of 

(11. 5) ,  the  (equally  appropriate)  procedure  of  regarding  y  to  be  the 

s 

basic  parameter  —  given  in  terms  of  physical  quantities  via  (II. 5)  and 
the  text  Eq.  (2.57b)  will  be  followed.46 

The  next  step  is  to  develop  a  Taylor  expansion  of  F(y)  about  y 

s 

(up  to  and  including  terms  quadratic  in  y-y  ).  For  this,  F"(y  )  is  re- 

s  s 

quired.  Differentiating  (II. 4),  and  using  (II. 5)  (together  with  algebraic 
manipulations  of  the  type  employed  in  establishing  the  final  form  of 

(11. 5) ),  one  obtains 


F"(ys)  -f 


yo  y  +1 
s  s 


1/2,  .1/2 

ys 


(II. 6) 


Using  the  standard  steepest-descents  approximation  of  replacing 
F(y)  by  its  Taylor-evpansion  in  the  vicinity  of  the  saddle  point,  and 
extending  the  lower  limit  of  integration  to  minus  infinity,  one  has 
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where  F"(y  )  is  given  by  (II. 6),  and  where 

S 


'  'fv  m 
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(a)  an  attractive  polarization  type  interaction  (^1/r  ),  arising  from 
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22.  Cf.  Ryshik  and  Gradstein,  Tables  (Deutsche  Verlag  der  Wissenshaf ten, 
1963),  formulae  7.322  and  7.335. 


23.  In  Drude's  theory,  absorption  (in  the  wing-region,  wx  >>  1)  is  pro¬ 
portional  to  the  real  part  of  the  conductivity.  This  quantity  is  given 
by  the  formula 


a(w) 


e2  1 

2  i 


Mio 


(per  charge  carrier).  Corresponding  to  this  real  part,  there  exists 
a  nonvanishing  energy  transfer  from  the  electromagnetic  field  to  the 
charged  particle,  given  by 


Dividing  by  the  incident  flux,  c£2/2tt  one  then  has  for  the  absorption 
cross  section, 

2 

q  =  - —  ,  q.e.d. 

ttMcco  t 


With  regard  to  the  factor  4/3,  occurring  in  (2.52)  it  should  be 
pointed  out  that  the  factor  1/x  in  Drude's  expression  actually  denotes 
the  momentum-transfer  collision  rate,  and  should  thus  be  equated  to 
the  total  collision  rate  times  the  factor^l  -  cos  (where  is 
the  scattering  angle) .  Now  it  is  a  relatively  straight-forward  exercise 
to  show  that,  in  specular  reflection  from  a  plane  surface,  an  average 
over  all  possible  angles  of  incidence,  ^  cos  -  1/3;  hence 

^1  -  cos  ip  ^  =  4/3. 
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24.  The  circumstance  that  (2.55)  is  an  impact-type  expression  should  not 
be  surprising,  since  the  physical  significance  of  its  condition  for 
validity  (namely  a y  <<:  1)  is  simply  that  an  appropriately 

defined  "collision  time"  (in  this  case,  the  interaction  range,  1/a, 
divided  by  some  average  velocity,  (v^+v^)/2)  be  small  compared  to  the 
reciprocal  of  the  external  frequency,  u.  This  condition  is  one  of  the 
standard  criteria  for  the  applicability  of  impact-type  theories  (cf. 

T.  Holstein,  Phys.  Rev.  79 .  744  (1950). 

25.  For  this  table  kT  and  M  were  set  equal  to  300°K  and  twice  the 
hydrogen-atom  mass,  respectively. 

26.  As  remarked  in  the  introduction,  even  this  upper  limit  will  be  found 

to  be  substantially  smaller  than  the  experimental  value  for  H„0  far- 

3  1 

wing  absorption  reported  by  Burch. 

27.  This  formula  actually  gives  the  r.m.s.  relative  velocity,  but  is 
sufficiently  accurate  for  our  purposes. 

28.  Actually,  the  quoted  paragraph  gives  this  condition  as  (u)-<jjq)  tc  <  1, 
where  ojq  is  the  unperturbed  line  frequency.  In  our  problem,  as 
discussed  lengthily  in  Appendix  I,  this  frequency  (i.e.,  any  one  of 
the  dominant  rotational  frequencies)  is  set  equal  to  zero. 

29.  In  the  case  at  hand,  one  has  to  consider  the  perturbation  of  the 
rotational  frequencies  (or,  more  generally,  the  rotational  motion, 
itself)  arising  during  the  transient  existence  of  the  ^0-^0  collision- 
complex. 

30.  In  the  actual  physical  problem,  this  would  correspond  to  the  distance 
of  closest  approach  of  the  one  of  the  hydrogen  atoms  to  the  nearest 
atom  of  the  perturbing  molecule  (e.g.,  in  ^0-^0  collisions  the 
closest  distance  of  approach  of  one  hydrogen  atom  of  the  first  ^0 
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molecule  to  the  oxygen  atom  of  its  collision-partner).  Said  distance, 
as  well  as  the  magnitude  cf  V  ,  may  ultimately  be  estimated  by  various 
methods  (e.g.,  structure  and  hydrogen-bond  energies  in  ice  and  liquid 
water) . 

31.  This  statement  implies,  of  course,  that  one  doesn't  take  V  (x)  to  be 

a 

literally  constant;  otherwise  it  would  cancel  out  of  the  problem.  One 

has  always  to  understand  the  assumed  constancy  in  the  limiting  sense 

in  which  V  (»)  is  still  zero, 
a 

—2 

32.  The  characteristic  u>  dependence  of  (3.16)  of  classical  impact  theory 
manifestly  obtains  in  the  limit  "ft  to  -*  0. 

33.  Strictly  speaking,  we  should  set  V(rQ)  equal  to  the  initial  relative 
kinetic  energy  of  the  two  molecules.  However,  as  will  be  seen 
aposteriori;  the  individual  magnitudes  of  the  attractive  and  repulsive 
components  of  (3.20)  are  sufficiently  large  so  that  said  relative 
kinetic  energy  represents  but  a  small  correction. 

34.  Only  the  barest  outlines  can  be  given  at  this  time,  since  no  specific 
calculations  have  as  yet  been  carried  out.  It  is  hoped  to  have  something 
more  definite  for  the  final  report. 

35.  in  particular,  there  are  actually  two  potential  minima,  each  of  which 

are  characterised  by  the  formation  of  the  well-known  hydrogen-level 
configuration:  0  —  H  -  0. 

36.  cf.  G.  C.  Plmental  and  A.  L.  McClellan,  The  Hydrogen  Bond.  (Freeman 
and  Co.,  1960),  p.  126. 

37.  The  justification  for  the  use  of  this  theory,  as  well  as  all  other 
developments  of  the  hindered-rotation  approach,  will  have  to  bs 
reserved  for  the  final  report. 
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It  is  herewith  once  and  for  all  assumed  that  the  (field-free) 
Hamiltonian  of  the  system  is  such  that  there  is  no  preferred  spatial 
direction.  Thus,  any  cartesian  component  may  be  used;  in  fact 
l<nl^zm>|2  will  eventually  be  replaced  by  -j  | <n| p |m> | ^. 

Literally,  s  arises  from  the  recipe  that  the  field  be  slowly  "turned 
on  (as  e  ) ;  as  is  known,  such  a  recipe  leads  to  a  correct  description 
of  dissipative  effects.  In  the  standard  Kramers-Heisenberg  formula, 
the  fact  that  the  energy  denominators  E  -  E  +  tiu)  are  nonzero  permits 
one  to  go  to  the  limit  s  =  0. 

(I. 5a)  is  of  course  violated  for  the  problem  of  interest  to  us,  but 
should  be  adequate  for  the  restricted  goal  of  this  appendix,  namely, 
a  study  of  the  adequacy  of  replacing  rotational  motion  by  linear 
motion. 

As  has  Keen  pointed  out  sometime  ago  by  E.  P.  Gross,  J.  Phys.  Chem. 

23.  p.  141.5  (1955),  this  latest  feature  of  "conservation  of  orientation" 
is  contradicted  by  the  "strong-collision"  model  of  van  Vleck  and 
Weisskopf,  Rev.  Mod.  Phys.  1948,  1946,  which  assumes  complete 
disorientaticn  per  collision.  Incidentally,  although  the  present 
treatment  differs  from  the  Boltzmann  equation  approach  of  Gross,  it 
leads  to  equivalent  results. 

In  a  stochastic  theory,  the  symbol<^. . . , acquires  the  meaning  of 
a  mixed  stochastic-thermodynamic  average,  the  stochastic  aspect 
consisting  of  averaging  over  all  stochastic  possibilities  with 
appropriate  weights. 

Actually,  in  the  case  of  water  vapor,  when  rotational  motion  is  that 
of  an  asymmetric  top,  three  rotational  constants  (corresponding  to  the 
three  moments  of  inertia)  exist.  The  value  given  here  is  that  which 
appears  to  be  the  most  relevant  for  rotation  of  the  dipole-moment. 
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In  these  equations,  the  subscript  "red"  is  dropped  for  notational 
convenience. 

Inspection  shows  that  the  r.h.s.  of  (II. 5)  is  monotonic;  there  is 
hence  only  one  value  for  the  saddle  point,  y  . 

The  value  yg  for  a  particular  value  of  the  parameter,  a,  may  ulti¬ 
mately  be  determined  by  any  convenient  numerical  or  graphical 
procedure. 


Chapter  2 


STIMULATED  RAMAN  EFFECT 


-*  > 

•*  i 


1. 


INTRODUCTION 


In  a  previous  report1  the  general  time  dependence  of  Raman  back- 
scattered  waves  was  considered  and  applied  to  the  specific  problem  of 
the  self-quenching  of  a  very  intense  laser  beam  in  the  atmosphere.  In 
this  report  we  shall  first  review  some  of  the  general  results  obtained 
previously.  Next,  we  shall  describe  the  action  of  a  Raman  amplifier. 
Finally,  an  attempt  to  quantitatively  describe  a  recent  experiment^  in 
Raman  pulse  generation  will  be  presented.  The  research  for  this  last 
description  has  not,  as  yet,  been  completed. 

Let  us  consider  the  interaction  of  a  beam  of  laser  light  passing 
through  a  length  of  Raman  active  substance,  with  a  beam  of  light  at  the 
first  Stokes  frequency  of  the  laser  light  in  that  substance.  If  the 
second  beam  is  traveling  in  a  direction  which  is  opposite  to  the  first 
it  will  be  amplified  by  the  stimulated  emission  of  backscattered  light. 
Assuming  that  all  dissipative  processes  -  including  the  stimulated 
emission  of  forward  scattered  Raman  photons-  are  inconsequential  (i.e., 
they  contribute  a  negligible  amount  to  the  depletion  of  the  laser  beam 
intensity  within  the  length  of  Raman  active  material  considered)  then 
the  densities  of  the  laser  photons  p£(x,t)  and  amplified  backscattered 
Raman  photons  p^(x,t)  in  the  substance  are  related  by  the  rate 
equations1 

(aT+c^r)pji  "  -  p*(apb+6)  (i) 

( aT  "  c  pb  =  +  pi(apb+6)  (2) 

where 

c3  8, 

a  =  o — cPT°(n)  and  6  =  cp  o(tt)  ~  . 

2v  Av  1  ^ 


*  % 

•<> 
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The  quantity  o(it)  is  the  Raman  backscattering  cross  section,  is  the  numb 
density  of  the  medium,  Av  is  the  line  width  of  the  amplified  wave  and 
is  the  effective  solid  angle  for  Raman  amplification  —  will  be 
discussed  in  greater  detail  later  in  this  report.  The  change  of  variables 


s  = 


T 


9 


R^(s,i)  =  P  ^  (x ,  t )  ,  F_k  ( s ,  t  )  =  pb(x,t)  , 


enables  us  to  write  Eqs .  (1)  and  (2)  as 


9R 

2  IT  ‘  -V°V6)  (3) 


9R. 

2~  "  R«>VS)  (4) 


The  general  solution  of  these  equations  is  given  by 


R£(s,t)  = 


IhL 


G(t)  -  f  F(s) 


(5) 


Rb(s,t)  = 


a  a  .  a 


G£H. 


(6) 


G(t )  -  f  F(s) 


•  • 

where  F  and  G  are  arbitrary  functions  and  F  and  G  are  the  derivatives  of 
these  functions  with  respect  to  their  arguments. 


Equations  (5)  and  (6)  have  been  obtained  previously.  The  appli¬ 
cation  of  boundary  conditions  to  these  equations  in  such  a  way  as  to 
describe  first  a  Raman  amplifier,  and  then  a  Raman  pulse  generator,  is 
the  subject  of  this  report. 


*  *■ 


2. 


THE  RAMAN  AMPLIFIER 


*■ 


Consider  a  situation  in  which  a  Raman  active  material  occupies 
the  slab  0  <  x  <  L.  Let  laser  light  of  constant  intensity  be  incident 
on  this  slab  from  the  left  so  that 


P£(°.t)  =  PD  (7) 

for  all  time.  At  t  =  0  the  slab  is  irradiated  on  the  right  by  a  beam 
of  light  at  the  first  Stokes  frequency  of  the  laser  lirjht  in  the  sub¬ 
stance.  The  photon  density  of  the  Stokes  light  is  given  by  R(t) : 

R(t)  =  pb(L,t)  ,  R(t)  =  0  for  t  <  0.  (8) 

For  the  purposes  of  this  section  we  shall  assume  that  the  spontaneous 
creation  of  backscattered  photons  in  the  amplifier  is  negligible,  so  that 
we  may  neglect  6  in  Eq.  (6).  Defining  the  quantity 


n(t) 


£ 

2 


t 


dt'  R(t’) 

)  —00 


(9) 


we  may  employ  Eqs.  (5)  and  (6)  to  rewrite  the  boundary  conditions  (7) 
and  (8)  in  the  form: 


G(t)  -  |  F(t) 

G(t  +  -) 
c 

G(t  +  “)  -  f  F(t  - 


(10) 


(ID 


Equation  (11)  may  be  written  as 


G(t  +  ^)  -  n(t)  G(t  +  ^)  = 


-  f  n(t)  F(t  -  -) 
^  c 


(12) 


'l 

'P 
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Defining  the  quantity 


8  "  f  po  *  (13) 

solving  Eq.  (10)  for  G(t),  and  substituting  the  results  in  Eq.  (12)  we 
obtain  the  differential  difference  equation  for  F(t) : 

F(t  +  £)  +  [g  -  n  (t)  ]  F(t  +  h  -  gn  (t) F( t  +  -) 

c  C 

=  -  gn(t)F(t  -  ~)  (14) 

Ihis  equation  may  be  concisely  solved  in  terms  of  the  quantity 

E(t)  =  F(t)  egt  .  n5) 


Equations  (5)  and  (6)  imply  that  we  may  arbitrarily  choose  E(0)  equal  to 
unity.  Equation  (14)  may  be  rewritten  in  the  form 


E  (t) 


1  +  6  (t  -  i) 

c 


dt'Vt")  e8t"  + 


r 

1-8 


dt  *E(t *-  i)e"gt’  -n(t,) 

_  C 


(16) 


where  0(t)  ■  1  for  t  >  0  and  6(t)  -  0  for  t  <  0.  Equation  (16)  relates 
the  value  of  the  function  E(t)  at  a  time  t  to  its  values  at  times  less 

than  t  -  2  .  Since 

c 


E(t)  «=  1  for  t  i  — 

c 

we  may  step-wise  extend  our  knowledge  of  E(t)  in  time  intervals  of 
duration  2  -  by  first  inserting  Eq.  (17)  into  Eq.  (16)  to  give 


... 
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9  L 

r  2  87 

E(t)  =  1  +  0(t-  -)  e  C 

c 


n(t"  7>  + 


dt"  n(t")ft  dt'n(t')  eg(t"-t’)+n(t")-n(t') 

J  rt 


for  t  <  3  —  ,  (18) 


and  then  inserting  Eq .  (18)  into  Eq.  (16)  to  give  E(t)  for  t  <  5  -  etc 

~  c  * 

The  laser  and  Raman  photon  densities  at  any  point  in  the  amplifier 
are  related  to  E(t)  in  the  following  manner 


<4k 


Rl(s-T>  -  po 


!(s) 


E(s)  + 


df  e8(s-t ')  i(tl) 


(19) 


Rb(s,  x) 


2  -28T 


E(x) 


E(s)  + 


dt’  e^8'^  E(t') 


(20) 


The  amplified  signal  is  given  by 


Pb(o,t) 


2  E(t) 
a  E(t)  * 


(21) 


f 


Equations  (16)  through  (21)  describe  the  behavior  of  the  amplifier 
for  an  arbitrary  pulse.  As  an  illustrative  example  let  us  consider  an 
input  square  pulse  of  infinite  duration  and  photon  density  pi  i.e.,  let 

2 


pb(L,t)  "  7^(t)  ■  fh0(t)  -  p1e(t), 

where  h  -  ap1/2  by  definition. 


(22) 
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If' 


In  this  case  Eq.  (18)  gives 


E(t)  =  1  + 


(g+h)(t-  p 


(g+h)‘ 


-l]  +  (g+h)g(t-  ■£■) 


for  L/c  <  t  <  3  L/c,  and  Eq.  (21)  yields 


Pb(0,t)  =  pl 


(g+h) 


(g+h)  (t-  f, 


■) 

+  g. 


h{h[f(8Fh' (t~  c)-0+  (g+h)(t-  ^)J  +  (g+h) 


2  p-2el 


in  that  same  interval.  Equation  (24)  exhibits  the  gain  and  "pulse-width 
sharpening"  characteristic  of  laser  amplifiers.  Expanding  the  amplified 
signal  about  its  time  of  ai rival  (L/c)  we  find  that 

pb(°’t+c)7~M^  PL  G{1  -  (G-l)ht  +  ...}  ,  (25) 


where 


L 

a—  p 
c  o 


is  the  amplifier  gain. 

The  step-wise  integration  of  Eq.  (14)  indicated  by  Eq.  (16)  is 
useful  for  pulses  which  are  no  longer  than  several  times  2  L/c  (note, 
for  L  =  30  cm. ,  2  L/c  =  2  nanoseconds) .  For  very  long  pulses  the  large 
number  of  integrations  required  to  describe  the  amplified  pulse  at  late 
times  would,  of  necessity,  produce  considerable  difficulties.  However, 
for  the  input  square  pulse  of  infinite  duration  described  by  Eq.  (22) 
some  simplifications  can  be  made  and  an  analytic  expression  for  E(t) 
for  all  times  may  be  obtained.  Substituting  Eqs.  (15)  and  (22)  into 
Eq.  (14)  we  obtain 

2g- 

E(t  +  ~)  -  (g+h)E(t +“)  =  -gh  e  c  E(t-  £)  (27) 
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with 


E(t)  =  1  for  t  <  and  E(^) 


h  e 


(28) 


Defining  the  Laplace  transform  of  E(t+  £■)  by 


E(p)  = 


dt  e“ pt  E(t+  £) 


(29) 


we  readily  obtain 


,  2gi 

E(p)  =  i  +  he  C 

P  P  D(p  i 


(30) 


where 


2  2~(g-p) 

D(p)  =  p  -  (g+h)  p  +  gh  e 


(31) 


We  may  use  the  Bromwich  integral 


E(t+  b 


d+i°° 

if  pt  • 

2ttT  j  dP  e  E(P>  . 
d-i°° 


(32) 


where  d  is  such  that  the  contour  lies  to  the  right  of  all  the  poles  of 
E(p),  to  express  E(t  +  £)  in  terms  of  the  finite  sum 

2g~  ^  2—  n 

E(t+  ~)  -  1  +  h  e  ^(-gh)n  e  An(t-2^-  n)  ,  (33) 

n«*0 


ct 


f°r  N<2^-<N+1,  N  *  integer, 


where 


,n+l 


An(t) 


(P-g)  Pt 


(n+1) !  n+1  .  n+1 

dp  lp-g-h) 
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] 


p=o 


4  f  ^  ePt  1 

n:  dPn  LPn+2  J 


P=g+h 


(33) 


(34) 


This  result  is  completely  equivalent  to  that  obtained  from  Eq.  (16). 

It  is  much  simpler  to  obtain  and  express  analytically,  but  still  quite 
difficult  to  employ  for  N  >>  1.  In  that  case  E(t)  can  be  obtained  from 
Eq.  (32)  by  appropriately  distorting  the  integration  path.  For 
2  L/c  i  tQ  =  h  "'"-g  the  result  is 


E(t  +  -)  - h 

c'  t  00 


2&Z  P0“g  P0t 

0  - ■ -  0 

P0D  (P0) 


(35) 


where  pQ  is  the  pole  of  (p-g)/D(p)  with  the  largest  real  part.  It  can  be 
shown  that  pQ  is  on  the  real  axis  and  is  a  zero  of  D(p)(but  not  the  zero  p=g) . 
The  pole  pQ  is  obtained  graphically  (see  Figs.  Ha  and  lib)  by  plotting 
the  parabola  (g+h)p-p  and  the  exponential  gh  exp[2L(g-p)/c]  versus 
p  for  p  real.  These  curves  intersect  twice:  once  at  p=g  and  once  at  p=p  . 
Examining  these  curves  we  find  that  for  2L/c  <  tQ,  0  <  p  <  g,  and  for 
2L/c  >  tQ,  g  <  pQ<  g  +  h.  Using  Eq.  (35)  in  Eqs.  (19)  and  (20)  we  obtain 


p  (x,t)  - - »  p 

J6  t  ->  oo  O 


g-P, 


g-P0  e 


“2(g~po)x/c  ’ 


(36) 


p,  (x,t)  - >  j-  p 

D  ^  oc  &  rO 


p0-g 


P0-g  e 


2(g-PQ)x/c 


(37) 


from  which  with  the  aid  of  D(p^)  =0  —  it  is  easy  to  show  that 


(38) 


(,b(0’t)  rn?  |p„  •  09) 

Equations  (36)  and  (37)  -  describing  p£  and  pfe  in  the  interior  of  the 
amplifier  after  a  very  long  time  -  are  graphically  displayed  in  Fig.  III. 
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Equation  (39)  indicates  that  although  the  leading  edge  of  the  input  pulse 
described  by  Eq.  (22)  may  be  amplified  to  a  value  greatly  in  excess  of 
Pq+  the  tail  of  the  pulse  is  always  given  by  2pQ/a  <  pq+  p^.  The 

sum  of  Eqs .  (38)  and  (39)  leads  to  the  statement  of  conservation  of 
energy:  p£(L,t)  +  pb(0,t)  =  pq+  p^ 
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3. 


THE  RAMAN  PULSE  GENERATOR 


In  this  section  we  shall  employ  a  model,  based  on  Eqs.  (5)  and  (6), 

Q  f 

and  the  theory  of  "self-trapping,"  *  to  quantitatively  describe  the 

r\ 

results  of  an  experiment  performed  by  Maier,  et  al.  In  this  experiment 
a  carbon  disulfide  cell  was  illuminated  by  a  beam  of  laser  light.  Ob¬ 
servations  at  the  entrance  of  the  cell  revealed  intense  short  pulses  of 
backward  scattered  Stokes  radiation.  These  pulses  had  a  duration 
^  3x10  sec  and  a  peak  power  one  order  of  magnitude  higher  than  the 
incident  laser  power.  Since  the  gain  in  the  CS2  cell  was  quite  large, 
it  is  not  surprising  that  the  emitted  pulse  was  very  intense.  However, 
the  mechanism  for  the  production  of  the  initiating  pulse  of  backscattered 
Stokes  radiation  remains  to  be  explained.  The  obvious  candidate  for  such 
a  mechanism  is  spontaneous  emission.  However,  it  has  been  previously 
demonstrated  that  if  a  backscattered  Raman  wave  in  a  uniform  laser  beam 
was  generated  solely  by  the  spontaneous  emission  in  that  beam,  then  the 
intensity  of  the  Raman  beam  cannot  exceed  that  of  the  input  laser  beam 
(a  physical  explanation  of  this  phenomena  will  be  given  later).  Maier, 
et  al. ,  suggest  that  "A  probable  mechanism  for  the  initiation  of  the 
pulse  is  the  abrupt  onset  of  backward  stimulated  Stokes  emission  near 
the  exit  cell  surface,  accompanying  the  occurrence  of  laser  self-focusing 
in  that  region."  In  support  of  this  suggestion  they  report  that  in 
measurements  of  the  dependence  on  cell  length  of  the  laser  threshold 
power  for  pulse  formation,  they  found  a  result  which  was  characteristic 
of  the  self-focusing  effect.  As  further  support  for  this  suggestion 
we  shall  construct  a  specific  model  of  the  influence  of  self-focusing 
on  the  amplification  of  spontaneously  emitted  Raman  waves.  The  validity 
of  our  model  may  be  measured  by  the  degree  of  success  we  attain  in 
predicting  the  height  and  shape  of  the  emitted  pulse. 

Consider  first  a  Raman  active  medium  and  a  beam  of  laser  light  in¬ 
cident  on  it  from  the  left.  Let  the  photon  density  of  the  incident  beam 


72 


<+ 

Jr 


be  the  constant  p^  and  let  the  time  of  incidence  be  t  =  0.  The  boundary 
conditions  are 


P0(O,t)  =  p  0(t),  p,  (x,t)  =  0  for  t  <  0. 


Equation  (4)  indicates  that  is  a  continuous  function  of  s  and  hence, 
since  p^  =  0  for  t  <  0  we  have 

1^(0, t)  =  0  . 

Using  Eqs.  (40)  and  (41)  in  (5)  and  (6)  we  have"*" 


1  -(g+3/2)s 

p  (x,t)  =  p  e(s)8 — - -  ■■■  ~  ,e - r--...  .  (I 

b  °  (g+6/2)  e6x/<:  -g+g  e-(8+6/2)s 

Defining  tQ  to  be  the  time  required  for  the  medium  to  generate  a  Raman 
wave  at  the  entry  plane  which  is  half  the  intensity  of  the  incident  wave, 
i.e.  , 


Pb(0,to) 


2  Po  ’ 


we  obtain  the  relation 


9 

t  (o  )  =  — £n(— — ) 
o  o  ap  3 

o 

to  an  excellent  approximation.  From  Eq.  (44)  it  is  clear  that  t  de- 

o 

creases  with  increasing  p^;  ie.,  as  the  intensity  of  the  incident  wave 
is  increased,  the  time  required  to  generate  a  Raman  wave  equal  to  half 
the  intensity  of  the  incident  wave  is  decreased. 


Consider  now  an  arrangement  in  which  a  laser  beam  is  incident  from 
the  left  on  a  Raman  active  medium  which  occupies  the  space  x  >  0,  and 
propagates  through  the  medium  for  a  distance  L  maintaining  a  constant 
cross  sectional  area  (see  Fig.  IV).  At  x  =  L  the  beam  cross  section  is 
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assumed  to  decrease  radically  due  to  self  trapping.  Let  the  laser  and 
Raman  photon  densities  be  denoted  by  p  ^  and  for  0  <  x  <  L,  and 
p^2  and  P2  f°r  x  >  L.  For  conservation  of  flux 

p£2(L,t)  =  X  pu(L,t)  (45) 

by  hypothesis.  In  Eq.  (45),  X  >>  1  is  the  ratio  of  the  beam  area  in 
region  1  to  that  in  region  2.  Self  focusing  occurs  because  the  laser 
beam  is  sufficiently  intense  to  induce  significant  changes  in  the  index 
of  refraction.  These  changes  are  large  enough  to  induce  a  lens-like 
condition  in  the  region  about  x  =  L.  The  initial  laser  beam  in  region  1 
is  consequently  focused  into  a  beam  of  a  much  smaller  diameter  in  region  2. 
In  region  2  the  focusing  effect  of  the  increased  index  of  refraction  is 
exactly  compensated  for  by  the  defocusing  effect  of  diffraction  and  the 
beam  propagates  with  a  constant  cross  sectional  area.  Assuming  that 
dispersion  is  sufficiently  small,  we  may  apply  the  principle  of  optical 
reversibility  to  those  Raman  photons  which  are  spontaneously  emitted 
within  the  cone  angle  for  amplification  in  region  2.  That  is,  upon  in¬ 
cidence  from  the  right  on  the  region  about  x  =  L  they  will  be  defocused 
by  the  same  factor  as  that  for  focusing  the  laser  beam  incident  upon 
this  region  from  the  left.  These  considerations  lead  to  the  boundary 
condition 

P2(L,t)  =  X  p1(L,t)  .  (46) 

In  a  beam  of  uniform  cross  section  the  spontaneous  emission  is  produced 
at  a  uniform  rate  along  the  beam.  Consider  two  photons  spontaneously 
emitted  at  x^  and  x2  with  x^  <  x2»  As  they  travel  to  the  left  they  will 
be  amplified.  When  the  intensity  of  the  wave  emitted  at  x2  is  comparable 
to  the  initial  laser  intensity  the  wave  emitted  at  x^  also  has  an  in¬ 
tensity  comparable  to  the  initial  laser  intensity.  The  laser  intensity 
at  location  of  the  x,-wave  is  therefore  depleted  and  the  gain  is 

t. 

decreased.  Consequently,  a  saturation  effect  occurs  and  the  x2~wave  is 
prohibited  from  exceeding  the  initial  laser  wave  in  intensity.  Since 
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O  x^  and  x2  were  chosen  arbitrarily  the  above  argument  holds  throughout  the 

beam  and  the  backscattered  spontaneous  Raman  beam  cannot  exceed  the  initial 
laser  beam  in  intensity.  This  consequence  of  Eqs .  (1)  and  (2)  was  previously 
demonstrated  in  a  purely  mathematical  fashion.  In  a  beam  of  nonuniform 
cross  section  which,  for  example,  corresponds  to  the  self-trapping  ar¬ 
rangement  discussed  above,  the  backscattered  wave  in  the  region  x  >  L  will 
grow  to  the  laser  intensity  in  that  region  before  a  large  decrease  of  the 
laser  intensity  in  the  region  0  <  x  <  L  can  occur.  This  fact  is  a  clear 
consequence  of  Eq.  (44)  and  the  fact  that  the  laser  intensity  is  much 
greater  in  region  2  than  region  1.  The  backscattered  wave  will  propagate 
through  the  plane  at  x  =  L  and  will  reappear  to  the  left  of  this  plane 
with  an  intensity  comparable  to  the  laser  intensity  there.  Since  the  gain 
in  region  1  is  as  yet  undiminished  the  backscattered  wave  will  now  be 
amplified  to  an  intensity  far  in  excess  of  the  initial  laser  intensity. 

The  physical  argument  given  above  is  readily  translated  into  a 
quantitative  description  of  the  Raman  pulse.  Eqs.  (3)  and  (4)  are  valid 
in  regions  1  and  2  separately  if  they  are  modified  to  include  the  fact 
that  the  constant  "0"  is  different  in  the  two  regions.  Consequently, 

R  and  R^  may  be  expressed  in  terms  of  F^,  and  0^  through  Eqs.  (5) 
and  (6).  This  is  also  true  for  R£2,  R2>  F2,  G2  and  0^  The  constants 
"0"  differ  by  the  factor  ft  -  the  effective  solid  angle  for  Raman  ampli¬ 
fication.  The  effective  solid  angle  is  obtained  in  the  usual 

manner'*'  for  a  baam  of  length  L  and  a  given  cross  section.  The  effective 
solid  angle  "fl2"  corresponds  to  that  of  a  cone  with  a  half  angle  equal 
to  the  angle  of  total  internal  reflection  in  region  2.  Since  there  are 
four  independent  functions  (F^  2  and  G^  2)  four  boundary  conditions  are 
required.  Two  boundary  conditions  are  given  by  Eqs.  (45)  and  (46). 

A  third  condition  describes  the  source  function.  For  simplicity  we 
may  take 

pu(°,t)  =  pQ  6  (t )  .  (47) 

n 
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The  final  condition  is  that  there  is  no  radiation  at  the  frequency  of 
the  Raman  waves  incident  on  the  system  from  the  right;  i.e., 

R2(0,t)  =  0  .  (48) 

Using  Eqs .  (45)  through  (48)  in  the  modified  Eqs.  (5)  and  (6)  leads  to  an 

expression  for  p^(0,t)  in  terms  of  a  function  which  is  the  solution  of  a 

differential  difference  equation  as  in  the  preceding  section.  Since  the 

4 

theory  of  self-focusing  gives  us  a  relationship  between  p  and  L,  it  is 
possible  to  determine  the  height  and  shape  of  the  exciting  pulse  solely 
in  terms  of  the  incident  laser  intensity  and  X. 

The  program  outlined  above  for  describing  the  backscattered  pulse 
has  not  as  yet  been  carried  out  because  of  some  difficulties  in  solving 
the  differential  difference  equation.  As  an  aid  in  gaining  insight  into 
the  nature  of  the  solution  to  the  above  boundary  value  problem  a  simpler 
problem  has  been  solved.  This  simpler  problem  is  a  mixed  boundary- 
initial  value  problem.  For  this  latter  problem  all  space  (i.e., 

-  «,  <  x  <  oo)  is  filled  with  Raman  active  material.  The  seK  focusing 
occurs  at  x  =  0  (in  contrast  to  x  =  L  in  the  preceding  problem)  so  that 

P£2(0,t)  =  X  p£i (° . t )  ,  p2(0,t)  =  X  P;L(0,t)  .  (49) 

Initially  the  Raman  and  laser  waves  are  arranged  so  that 

pU(x’0)  =  p0  »  Pjt2^X,°^  =  pi(x’°>  =  P2(x»°)  =  0  •  (50) 

Equation  (48)  is  valid  in  this  problem.  An  additional  condition  is 

R^l(s,0)  continuous  in  s  .  (51) 
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To  solve  the  modified  Eqs .  (5)  and  (6)  we  first  define  the  quantity 


(ap  +  B,)L/c 

=  e  and  x  =  t  -  L/c  >  0.  Examinations  of  Eq.  (53)  both 

analytically  and  with  the  aid  of  a  computer  show  that  iKt)  drops  very 
sharply  -  like  exp(-Xgt)  -  from  its  initial  value  until  it  reaches  a 
minimum  at  t1  =  30  (Xg)-1.  It  increases  gradually  from  that  point  onward. 

The  backscattered  photon  density  p.  exhibits  the  expected  pulse-like 

-i 

character.  It  rises  to  a  maximum  at  a  few  times  (Xg)  earlier  than  t^. 
The  maximum  ranged  between  3p^  and  15CpQ  for  the  parameters  considered, 
namely,  10  cm  <  L  <  100  cm  and  10  <  X  <  150. 

It  is  our  intention  to  complete  the  analysis  of  the  pure  boundary- 
value  problem  during  the  remainder  of  this  contract. 
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Figure  1.  A  Raman  amplifier  of  length  L  laser  light  — 
photon  density  pQ  —  is  incident  from  the  left 
at  x  =  0.  The  signal  —  light  at  the  first  Stokes 
frequency  with  a  photon  density  p^(L,t)  -  is 
incident  from  the  right  at  x  *  L.  The  amplified 
pulse  —  photon  density  p^(0,t)  —  is  emitted 
at  x  ■  0, 
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Figure  3.  The  photon  density  within  the  amplifier 
in  the  limit  t  ■+■  °°  . 
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Figure  4.  The  Raraar  pulse  generator.  A  Raman  active  medium 
occupies  the  space  x  >  0.  Laser  light  -  photon 
density  is  incident  on  this  medium  from  the 

right.  It  propagates  from  x  -  0  to  x  -  L  with  a 
constant  cross-sectional  area.  At  x  =  L  self 
focusing  occurs  and  the  area  of  the  beam  decreases 
by  a  factor  of  (1/A)  «  i.  The  laser  and 
photon  densities  are  and  p^.  respectively  for 
0  <  x  <  L,  and  p ^  and  p2  for  x  >  L. 
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formation  in  the  paragraph,  represented  as  (TS),  (S),  (C),  or  (V). 

There  le  no  limitation  on  the  length  of  the  abstract.  How¬ 
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l-  -^K*Y.WORPB:-  **7  wor<J*  technically  meaningful  terms 
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fiera,  such  aa  equipment  model  designation,  trade  name,  miiitarr 
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